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IMP: INTERACTIVE MASS PROPERTIES PROGRAM 

Volume I - Program Description 
by 

William A. Stewart - Sigma Corporation 

SUMMARY 

This document describes a method of computing a weights and 
center of gravity analysis of a flight vehicle using interactive 
graphical capabilities of the Adage 340 computer. The equations 
used to calculate area, volume and mass properties are the 
equations derived in reference 1. The input/output methods are 
different than those presented in reference 1 because of the 
graphic support of the Adage computer and the geometry definition 
used for storing data. Several interactive program options are 
available for analyzing the mass properties of a vehicle. These 
options are explained in this report. An overall picture of the 
program with its displays and options is shown in figure 1. 

The document also describes the method and computer techniques 
for the calculation of the mass properties, c.g. location, 
enclosed volume, wetted area and planform area of any closed 
structure that has a plane of symmetry. The vehicle is described 
to the computer program by ordered sets of X, Y, Z coordinates 
of points on its surface. The X, Y Z coordinates are converted 
to quadrilateral elements for analysis. The mass properties of 
each quadrilateral may be computed from a thickness and density 
input for each quadrilateral or from a weight per unit area in- 
put at each point or from a combination of both. The weight 
per unit area can be composite of the real wall including skin, 
insulation, ribs, stringers, standoffs, brackets, etc. 

The mass properties of the quadrilateral elements are accumulated 
for eciCh section of the vehicle and for the total vehicle. The 
computed mass property totals may contain not only the contribu- 
tion from the distributed mass on the vehicle surface wall, but 
additional masses may be added as "black boxes" by specifying 
each one's c.g. location and mass properties. The added black 
boxes may be inside or outside the surface and do not have to 
be symmetrical with respect to X-Z plane. 
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FIGURE 1 INTERACTIVE MASS PROPERTIES PROGRAM 
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INTRODUCTION 


The estimation of the mass properties of a vehicle is one of 
the most important considerations in the design process and yet 
one of the most inexact engineering endeavors. While the calcu- 
lation of aerodynamic, propulsion and mission performance are 
based on widely recognized mathematical prediction techniques, 
the estimation of weight must be based largely on historical 
data. The art of weight estimation has evolved through the 
years by the diligent collection and correlation of component 
weights of previously built vehicles. New design weights are 
predicted on the basis of the component weights of past designs. 
Little information is usually available on the other properties 
such as volume, area, center of gravity and inertia of the com- 
ponents. Even for "weights only" calculation, historical methods 
lack validity in predicting advanced state-of-the-art design data. 
Further, the influence of small design perturbations eludes such 
gross prediction methods. 

During the later design phases when detail structured analyses 
are available for some components, the weights engineer must 
still revert to empirical analysis methods for subsystem and 
secondary structure weights such as wing leading edges, canopies, 
fairings, etc. 

Existing programs such as references 2 through 6 lean heavily 
on the historical approach for weight estimating. To obtain 
more reliable information, the designer must turn to structural 
analysis and subsystem design programs such as references 7 
through 13. The detailed geometry descriptions required for 
these classes of analysis are usually not known until the later 
design phases. At this later phase the design is usually frozen 
and little geometric perturbation is allowed. Further, small 
perturbation analysis using these programs is very expensive in 
terms of manpower and computer resources. 

Although the precise weight and other properties of the vehicle 
components are illusive, or at best a very complex calculation, 
mathematical equations for combining the components mass 
properties into total vehicle mass properties are quite concise. 
Mass property evaluations for the total vehicle as well as 
perturbations due to added or deleted masses are easily determined 
once the mass properties of all the components are known (or 
assigned) , 
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The significant contribution of the IMP program development to 
the general problem of computer aided design is threefold. 

First,, the program offers a unified approach to the accounting 
of the vehicle component mass properties. The fact that a com- 
ponent weight is not well known does not exclude it from being 
analyzed in a mass properties evaluation. Furthermore, detailed 
mass properties of some vehicle components may be determined 
elsev/here and still be included in the IMP analysis. Second, 
the distributed mass surface model incorporated in IMP allows 
a very convenient means of producing vehicle geometric perturba- 
tions early in the design process with some degree of confidence 
in the relative effect of the perturbation. As more detailed 
information on the design becomes available, the component input 
to IMP may be gradually replaced without distribution of the 
overall continuity of the weights analysis. Third, the IMP 
program allows the user to directly interact with the mass 
properties evaluation process through use of the Adage 340 com- 
puter. This highly sophisticated display computer allows the 
user an additional dimension in his analytical capabilities. 


VOLUME, AREA AND MASS PROPERTIES ANALYSIS 


The mathematical modeling of the IMP program permits the calcula 
tion of the mass properties, c.g. locations, enclosed volume, 
wetted area and planform area of any closed structure that has 
a plane of symmetry, c.g. fuselage, stiffened fuel tank, etc. 

The vehicle is described in a symmetrical manner and consists 
of sets of X, Y Z coordinates on its surface. The surface 
coordinates are converted to quadrilateral elements but since 
the four corners of the quadrilaterals are not necessarily 
coplanar, each quadrilateral is analyzed as two triangles. The 
mass properties are computed from a thickness and density input 
for each quadrilateral, from a weight per unit area input at 
each point, or from a combination of both. The weight per unit 
area can be a composite of the real wall including skin, insula- 
tion, ribs, stringers, standoffs, brackets, etc. 

The elemental mass properties are accuraualted for each section 
and for the total vehicle. The computed mass property totals 
containing the contribution from the distributed mass in the 
vehicle surface wall may be combined with additional masses 
specifying each one's c.g. location and mass properties about 
its c.g. The added black boxes may be inside or outside the 
surface and do not have to be symmetrical. 

The properties of the vehicle computed by the program are: 

1. Weight. 

2. Coordinates of the c.g. 

3. Three mass moments of inertia plus three products of 
inertias . 

4. Surface area. 

5. Enclosed volume. 

6. Projected area. 

These properties are output for each segment along with cumula- 
tive totals and complete vehicle values. 
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.Surface Model 


The numerical model of the surface shape to be analyzed is 
described by. quadrilateral elements. Each quadrilateral 
element consists of a grouping of four surface points. An 
organized set of quadrilateral elements form a segment of 
the configuration. A number of segments may be used to 
provide a complete description of the configuration. 

The configuration is usually (but not necessarily) positioned 
with its nose at the coordinate system origin and with the 
length of the body stretching in a positive X-direction. 

The coordinate system used for this program is the right 
hand Cartesian system as illustrated below: 


Z 



The configuration is assumed symmetrical about the X-Z 
plane, and the analysis accounts for the plane of symmetry. 

The configuration may consist of an arbitrary number of 
wing and body sections. The body sections are described 
in the positive half plane perpendicular to the plane of 
symmetry. The wings are described as airfoil sections 
parallel to the plane of symmetry. Canards and horizontal 
Btcibilizers are described similar to' wings. Fins cannot 
be described as airfoil sections but they may be described 
as body sections. 



Points along the X-axis are referred to as station coordinates 
A Y-Z plane is a station plane « The Y-Z coordinates in the 
station plane form a station contour * 



A portion of the mathematical surface model is presented 
in the illustration above. The X,Y,Z coordinates of the 
points where the solid grid lines cross are described to 
the program. Each point (I/J) is associated with a station 
number J and a contour point I, Thus, the circled point 


is 8,5 and it has coordinates 




The model 
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groups the surface points into quadrilateral elements which 
also have unique If J identification. The crosshatched one 
is 8,5 corresponding to its lowest numbered corner. The 
dashed diagonals illustrate the division of each quadri- 
lateral into two triangles for analysis. 

Airfoil sections are described in a similar manner as the 
body sections above. In the case of airfoil section, the 
parts are described in X.rZ planes parallel to the plane of 
symmetry . 


The objective in the analysis is to calculate the properties 
of each quadrilateral (two triangles) elements and trans- 
form the results back to the reference coordinate system, 
then accumulate the properties from all elements of the 
surface. 

Mass Properties of a Triangular Surface Element . - The 
triangular surface element illustrated below with its base 
a station J and apex at J+1 is assumed to be of uniform 
thickness and homogeneous. 

y 



As viewed from the outside of the vehicle, the local x,y,z 
coordinate system has its origin at (I,J) . The x-y plane 
is coincidental with the plane of the triangle. A is the 
height of the triangle. B is the base of the triangle. 



C is the offset of the apex from the y-axis. H is the 
thickness which is assumed to be small compared to A# B 
and C. 

The area, S of this elemental triangle is: 

S = AB/2 (1) 


The first moments of inertia for the same triangle are: 

( 2 ) 

(3) 


•S = PHS(B+C)/Z 
3C 


Sy = PHSA/2 
= 0 


(4) 


The second moments of inertia for the elemental triangle 


are: 


Ixx = pHA^S/6 

= pH(B^+BC+C^)S/6 


I =1+1 
zz XX yy 


(5) 

( 6 ) 

(7) 


Where p is the density of the thin triangle. Finally, the 
products of inertia are: 


I ^ = pH (B+2C) AS/12 
xy 

I =0 
yz 


"xz - ° 


( 8 ) 

(9) 

( 10 ) 


The coordinates of the centroid (center of gravity) of the 
triangular element are: 


X = (B+O/3 
y = A/3 


( 11 ) 

( 12 ) 
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The transformation of the inertias and products of inertia 
to the centroid yield: 



= SA^/18 

(13) 

^yy 

= X(B^-BC+C^)/18 

(14) 

^xy 

= SA(ZC-B)/36 

(15) 

T 

= 

(16) 

zz 

XX yy 


Transformation Matrix for Triangle Elements ^ - The mass 
properties above are given for a simple triangular element 
in the local y:,y,z coordinate system. The elemental mass 
properties must be transformed to the reference coordinate 
system X,Y,Z for cummulation with other elements. The 
transformation from the primary Y,.,Y,Z coordinate system to 
the local system is given by: 


1 

8 

= [dc] 

1 

0 O 
X X 

1 1 
X 

1 

L*J 

. 

.z - zoj 


(17) 


Where XO,YO,ZO are the coordinates of the prime system 
origin with respect to same arbitrary reference point and 
DC is a matrix of direction cosines. The elements of DC 
are: 



"XPX 

YPX 

_ZPX 


XPY 

YPY 

ZPY 


XPZ“ 

YPZ 

ZPZ. 


(18) 


For a unit vector in the x direction, XPX is the projection 
of X on X , XPY is the projector of x on Y , etc. The 
vertices of the triangular element in local x,y,z 
coordinates may be considered as vectices B and D, drawn 
from the origin (0,0,0). The components of B and D are 
developed as follows. That is^BX, is the X component of 
B and DY is the Y component of D, etc. Therefore: 



If: XO = X(J) (19) 

YO = Y(I,J) (20) 

ZO = Z (I, J) (21) 

Then: BX =*0.0 (B lies in YZ Plane) (22) 

BY = Y(I+1,J) - YO (23) 

BZ = Z (I+l, J) - ZO (24) 

DX = X(J+1) - XO (25) 

DY = Y(I+1,J+1) - YO (26) 

DZ = Z(I-1,J-1) - ZO (27) 


The cross product of B with D is a vector S2 in the z 
direction. Its magnitude is twice the surface area of 
the triangle. Let SX/SY and SZ be the components of S2. 
Since BX = 0, the cross product yields: 

SX = BY^DZ - DY^BZ (28) 

SY = BZ*DX (29) 

SZ = -DX*BY (30) 

Let B, D, S2 represent the magnitude of their respec^ve 
vectors. Then 'B'/B = unit vector in x direction and S2/S2 
= unit vector in z direction. 

It follows: 

S2/S2 cross ¥/B = unit vector in 


y direction (31) 

The components of each unit vector are a row in the matrix 
DC* Therefore: 

XPX = BX/B =0.0 (32) 

XPY = BY/B (33) 

XPZ = BZ/B (34) 

ZPX = SX/S2 (35) 

ZPY = SY/S2 (36) 

ZPZ = SZ/S2 (37) 

YPX = ZPX*XPZ - XPY*ZPZ « DX/A ' (38) 
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YPY = XPZ*ZPX 
YPZ = ZPX*XPY 


(39) 

(40) 


The preceding steps considered a triangle with its base on 
J and apex on J+1. The other half of each quadrilateral 
has its base on J+1 and vertex on J. The vertices are 
identified as follows: 

(I, J+1) 



The analysis for the inverted triangle is similar to the 
upright triangle. The only difference is in equations (13) 
through (21) . For the inverted triangle, 


XO = X(J+1) (41) 

YO = Y(I+1,J+1) (42) 

20 = Z(I+1,J+1) (43) 

BX = 0.0 (44) 

BY = Y(I,J+1) - YO (45) 

BZ = Z(I,J+1) - ZO (46) 

DX = X(J) - XO (47) 

DY = Y(I, J) - YO (48) 

DZ =* Z(I,J) - ZO (49) 


The above analysis is similar for airfoil type sections 
described in the x-z plane. 




Mass Properties of the Surface. - The c.g. and mass pro- 
perties expressions for a triangle are given in the 
previous section. The triangular element parameters A 
and B have been developed above and are always positive. 

2 2 2 

Now C =• D - A but C can be negative. However, C and 
its sign are given directly by the dot product of D with 
the unit vector in the x direction. 

C * DY*XPY+DZ*XPZ (50) 

The resulting c.g. coordinates are in the x, y, z system 
and the inertias are about axes through the c.g., parallel 
to the prime system. They are transformed to the system 
reference axes as follows: 

xo' 

YO + [DCl^ * 

ZO 




This is readily expanded since the inverse of DC is 
identical to its transpose. The products and moments of 
inertia parallel to the prime axes are assembled in PMIP: 



^XX 

”^XY 

-I T 

xz 

PMIP = 

^XY 

-lyY 

”^YX 


^ ^XZ 

■*^YZ 

"^zzJ 


(52) 


The minus signs on the products of inertia are a consequence 
of the coordinate system definition. This places the matrix 
PMIP in tensor form. It can be shown that the inertias are 
rotated parallel to the system axes by: 

[PMI]' = [DC]*^ * [PMIP]* [DC] (53) 

PMI contains the inertias about axes parallel to the system 
reference axes, but centered on the triangld s c.g. Trans- 
lation to the system origin makes use of the parallel axis 
theorem along with XBR, YBR, ZBR computed above, and is 
included with the simimation procedure described in the next 
paragraph. The products of inertia contained in PMI are 
the negative of their engineering definition, like they 
are in PMIP. 
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In the IMP analysis, all values in PMI are twice the 
true value since S2 = twice the area of the triangle, S2 
has been retained through all the computations. But each 
triangle has a symmetrical counterpart on the other side 
of Y = 0 whose contribution to the total is either equal 
or equal-and-opposite. In the latter case, the resultant 
is zero; i.e. , 1^^ = I^^ “ 0**^* 

Area Calculations . - The surface area and projected area 
result from the summation of the elemental areas from the 
£d?Ove cuialyses. The surface area is: 

A 1/2 E S2^ (54) 

S2 is twice the elemental area and the projected area (on 
the y-z plane) is; 

A = 1/2 Z SY. (55) 

P ^ ^ 

Volume Calculations . - The elemental volume of a surface 
element is the projection of the elemental area on the 
x-y plane multiplied by the average distance to the c.'g. 
of the element. The total volume is: 

V = E X . • ZPX. ‘ A. (56) 

i i 

The sign of the direction cosine accounts for whether 
the elemental volume is added or subtracted. 

The center of volime is simply: 

CV = I E ZPX^ 

The above calculations are approximate because of the 
assumption that the c.g. of the element is an adequate 
approximation for the entire elemental surface. The 
assumption is reasonable for small elements but may 
not be valid for a gross panelling of the configuration. 
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Black Box Contributions to the Mass Properties 

The surface model described above is suitable for describing 
shell structures such as tanks, skins, etc. However, the 
addition of mass not conforming the above model may be 
necessary in the complete description of the vehicle mass 
properties. In the IMP analysis, additional mass sources 
may be added by specifying each ones center of gravity (c. 
g.) location and mass properties. These mass sources called 
"black boxes" may lie inside or outside the surface and do 
not have to be symmetrical with respect to the x-z plane. 

The IMP program combines the detailed shape inputs and the 
black box data to produce the overall mass properties of 
the vehicle. 


Black box masses are described to the program in the 
following manner: 

1. Location of its c.g. in system reference 
coordinates. 

2. Total weight. ,• 

3. Moments and products of inertia about local 
axes through its c.g, 

4. Orientation of local axes relative to the 
system reference axes or relative to some 
intermediate, axes. 

The minimum required for a black box input is the specifica 
tion of items 1 and 2. If items 3 and 4 are omitted, 
the moments and products of inertia are assumed to be 
zero' and the orientation of the local axis is assumed 
to be coincidental with the system reference axes. The 
contribution of each black box is added to the cumulative 
totals produced by the analysis. Since the size and shape 
is not being specified, there is no addition to surface 
area or volume. Since the black boxes do not have to be 
symmetrical with respect to the x-z plane, 1^^^ 

0.0, and there can be a non-zero first moment with respect 
to Y. 

The moments and products of inertia of each black box are 
input directly to PMIP (see equation 52) . The program 
tadces care of filling the upper half and adding the minus 



signs. Rotation of the inertias uses equation 53. DC 
is obtained from input data as follows. Consider the 
local x,y plane: 


y 



Point £ is the origin. Points ]. and 2 are arbitrairy 
points on the +x and +y axes respectively. Coordinates 
of these three points^ measured in the system to which 
the inertias are being rotated, are inputs. It is a 
simple matter to construct unit vectors in the x and y 
directions and their components yield the first two rows 
in DC. equation 18. The third row is the cross product 
of these two unit vectors. 

In order to construct DC, it is not necessary that point 
0 be the c.g.. The coordinates x, y must be parallel to 
the axes to which the PMIP inputs are referenced, but 
•rotation is performed independently of translation. After 
the rotation of equation 53, the inertias in PMI are still 
centered on the c.g. of the black box. Only the rotation 
takes place and the contents of PMI are moved to PMIP. 
Thus, it is possible to make several rotations, if needed, 
to reach alignment with the system reference axes. The 
input to the program* can be developed from any convenient 
points on the black box and intermediate structure whose 
coordinates are known. The input does not always have 
to be in the same units as the system X,Y,5 data. 
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APPENDIX A 


SUBROUTINE DESCRIPTIONS 

This appendix describes the subroutines used by the interactive 
mass properties program and their general relationship to the , 
overall functions of the IMP program. Program MAIN is the main 
driver routine. Subroutine listings follow the descriptions 
below. 


• Subroutine ANRST 

ANRST is a Fortran subroutine which initializes variables used 
in the subroutines AN6 and AN7. This subroutine is internal to 
AN6 and AN7. 


Subroutine AN6 

AN6 is a Fortran subroutine that generates an input number with- 
in the input grid. 

Subroutine AN7 

AN7 is a Fortran subroutine that creates an input number within 
the input grid. This subroutine is used exclusively by WABCOM. 

The routine differs from AN6 in that AN6 is used by all other 
subroutines . 

Subroutine APNDF 

APNDF is an assembly language subroutine which allows the user 
to append data records from core to the end of a disk file, with- 
in a pack/volume. The subroutine is part of a general file 
access package called RANDOM 10. 

Image BKBX 

BKBX is an image subroutine written in Fortran which displays the 
mass properties of a black box and allows the user to make 
changes to the black box data interactively. This image is 
internal to subroutine WABCOM. 

Subroutine BLDMAN 

BLDMAN is a Fortran subroutine that converts an ATEXT file name 
to a format that can be used by the RANDOM 10 routines. 

Subroutine CGMOVE 

CGMOVE is a Fortran subroutine that supports the black box tracker 
option of the program. It allows the user to translate components 
with respect to a fixed fuselage. 
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Subroutine CLOSF 


CLOSF is an assembly language subroutine which supports the 
RANDOM 10 package. The subroutine is used to close a file which 
has been opened by the STRTF subroutine and was .filled by calls 
to APNDF. 

Subroutine CROSS 

CROSS is a Fortran subroutine that computes the cross product of 
two vectors . 

Subroutine DECFMT 

DECFMT is a Fortran subroutine that computes an output ATEXT 
format for an internal binary decimal number. The routine is 
used in conjunction with menu generation. 

Subroutine DSFIL 

DSFIL is an assembly language routine which is part of the 
RANDOM 10 package. The routine checks for validity or duplica- 
tion of a file name on a given pack/volume. If the file is 
found to already exist, the size and disk address are returned 
from the subroutine. 


Subroutine DSKER 

DSKER is an assembly language subroutine that supports the 
RANDOM 10 package. DSKER is an assembly language subroutine 
that provides a testing loop and restricts further operations 
until all requested disk operations are complete, and the disk 
and controller are ready for the next request. In general, the 
routine controls the timing within the RANDOM 10 package. 

Subroutine DSVOL 

DSVOL is an assembly language subroutine that checks the valid- 
ity and free space of a pack/volume. This subroutine supports 
the RANDOM 10 package. 


IMAGE GEO 

GEO is a Fortran IMAGE subroutine that displays the input values 
for the geometry calculations in the WABCOM subroutine. This 
subroutine allows the user to make changes to these inputs 
interactively. This image is internal to subroutine WABCOM. 
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IMAGE GRID 


GRID is an IMAGE that displays an alpha numeric grid that allows 
the user to make inputs interactively. 

Subroutine INPUT 

INPUT allows the user to construct a data file that is used 
interactively with the mass properties program. 

FUNCTION lOGT 

lOCT is a Fortran function that converts the decimal nximber to 
its octal equivalent. 


Subroutine LOAD 

LOAD is a Fortran subroutine that sets up array values so that 
they are properly used in the calculation of the WABCOM sub- 
routine. 


Subroutine LOCFLE 

LOCFLE is a Fortran subroutine that calculates the file sequence 
number within a pack/volume. This is a subroutine that takes a 
file name and locates its position within a pack/volume by its 
sequence number. 


Program MAIN 

This program is the executive control program for the IMP program. 

Subroutine MSM3X3 

This subroutine performs the multiplication of a three by three 
(3X3) matrix by another three by three (3X3) matrix. 

Subroutine MXD 

MXD is a Fortran subroutine that performs the multiplication of 
a -three by three (3X3) matrix by one by three (1X3) matrix. 

IMAGE' OPTIONS 

OPTIONS is a Fortran image subroutine that displays the various 
options used in the interactive mass properties program. This 
image subroutine allows the user to choose the option that he 
wishes. 
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Subroutine PAT 


PAT is a Fortran subroutine that performs the calculation of 
the parallel axis theorem. The parallel axis theorem is a 
method of summing center of gravities and inertias of several 
component parts to achieve the total values for the system. 

Subroutine RDSKOK 

RbSKOK is a Fortran language subroutine that verifies the 
existence of a pack/volume- and file on disk. If successful, 
the number of words and the beginning address are returned. 

Subroutine RESET 

RESET is a Fortran subroutine that initializes variables in the 
WABCOM subroutine. This subroutine is internal to WABCOM. 

Subroutine RESETl 

RESETl is a Fortran subroutine that initializes variables in the 
interactive mass properties program. This subroutine is internal 
to program MAIN. 

Subroutine RESET2 

RESET2 is a Fortran subroutine that initializes variables in the 
interactive mass properties program. This subroutine is internal 
to program MAIN. 

Subroutine RESETS 9 

RESET69 is a Fortran subroutine that initializes variables in 
the weights and balance subroutine. This subroutine is internal 
to WABCOM. 

Subroutine RNDRD 

RNDRD is an assembly language subroutine that randomly reads 
files access by the RANDOM 10 package. 

Subroutine SKETCH 

SKETCH is a Fortran subroutine that calculates the coordinates 
of the fuselage used in the black box tracker option. 

Subroutine STORE 

STORE is a Fortran subroutine that stores displayed data from 
the CRT onto a specified disk for later use. The subroutine 
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is called when the store option is selected from the screen by 
the light pen pick. 


Subroutine STRTF 

STRTF is an assembly language subroutine that creates or starts 
a new file in the RANDOM 10 subroutine package. 

Subroutine TOT 

TOT is a Fortran subroutine that sets up the accumulation array 
in the WABCOM subroutine. 


IMAGE TV7ERP 

TWERP is a Fortran IMAGE subroutine that displays the fuselage 
geometry generated by the subroutine SKETCH. 

Subroutine WABCOM 

WABCOM is a Fortran subroutine that calculates the mass properties 
of a flight vehicle with several component parts. There are two 
options in this subroutine, one being the geometry option which 
allows the user to input corner-point geometry and have the pro- 
gram output a black box component which is then summed into the 
total calculation of other black boxes. The other option is the 
black box option which allows a user to change the values of the 
mass properties of a component and reflect those changes to the 
total mass properties of the flight vehicle. 

Subroutine WABLST 

WABLST is a Fortran subroutine that generates the menus for the 
mass property data that is displayed on the CRT. 

IMAGE WABMNU 

WABMNU is a Fortran image subroutine that displays the numeric 
mass property data onto the CRT. 

Subroutine WABSET 

WABSET is a - Fortran subroutine that sets up data into mass 
property arrays that are used by the WABCOM subroutine. 

Subroutine XFORM 

XFORM is a Fortran subroutine that makes coordinate transforma- 
tions for the geometry option in WABCOM. 
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Subroutine XPOSE 


XPOSE is a Fortran subroutine that takes the transpose of a 
three by three (3X3) matrix, 

^ , Subroutine $NALTS 

$NALTS is an Adage system subroutine that halts the display of 
the present image. 

Subroutine $PPLCOPY 

$PPC.OPY -is an Adage system subroutine that makes a hard copy 
, of the present image displayed on the CRT. 

Subroutine $PTCH9 

This is an Adage system subroutine that specifies the CRT to be 
an input/output device. This subroutine needs to be called only 
once at the start of a program. 


Subroutine $ STOPS 

$STOPS is an Adage system subroutine that stops the execution 
of the program loaded and returns the control of the monitor 
system of the Adage. 
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REV 
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1 SUBROUTINE ANfa 

2 C Rl.O 

3 C CCNSTRUCTIOH OF ALPHA-NUMERIC INPUT DATA 

4 C 

5 COMMON-PARA- LFLG< 10 ) > MFLAGC4) . SUM . LETTER , NO . NFLG / INO/ NFLAG< 12) 

6 ♦ /JJJ>KKK/LLL.MMM.NNN,HAMEC2.60>.NME<2)/ IINUM.RRNUM 

7 . DIMENSION TSUM<20 > , SNC20 ) , ESUMC20 > . MASK<5) 

10 DATA TSUM-20>*=0 .-.SUM-0 LFLG<6)-0->LFLG<3>-l- 

11 DATA SN-20*0 .-.ESUM-20*0 .- 

12 DATA MASK-0 B. 77000000008. 7777000 000 B. 77777700008. 777777770 OB- 

13 IFCLFLGO) .EQ.O) GO TO 89 

14 100 CALL ANRST 
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VER 0 


REV 
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2 


AN6 


1 

2 

3 

4 

5 

6 
7 
10 
11 
12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 - 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

55 

56 

57 
60 
61 
62 

63 

64 

65 

66 


C A1 . 1 

C ONE DATA ENTRY PER LIGHT PEN DEPRESS 
C 

89 IF<LrLC<6) . EQ. 1> CO TO 89 
C A1.2 

C WHEN CLR IS PICKED BY LIGHT PEN 
C 

IFCLFLG<9) .EQ. 1) GO TO 100 
C A1 .3 

C CONSTRUCTION OF ALPHA WORD 
C 

IF<MFLAG<1 ) .EQ. 1 ) GO TO 60 
C A1 . 4 

C WHEN ENT IS PICKED BY LIGHT PEN 
C 

IF<LFLG<8) .EQ . 1 ) GO TO 999 
C A1 .5 

C DETERMINE IF NUMERICAL INPUT IS REAL OR INTEGER 
C 

;(f<:lflgc7>) 22 , 56.20 
20 ir<LFLG<5> .EQ. 1> GO TO 35 
K=K+1 

IFCLFLG<1 > . NE. 1 ) 60 TO 25 
22 ICNT=ICNT+1 
II=1-ICNT 
L=K+I I 

IFCICNT.EQ. 1> GO TO 24 
GO TO 30 

24 LFLG<7)=0 
GO TO 999 

25 II=K-1 

30 TSUM<K>=SUH 
SN<K) = 10 ,>»=*I I 
LFLG<7)=0 
GO TO 999 
35 M=M+1 

ESUM<M)=SUM 
LFLG<7>=0 
GO TO 999 
999 EXP=INUM=0 
RNUM=0 . 

IFCLFLG<5> . NE . n GO TO 43 

DO 42 1=1, M 

J=M-I 

EXP=EXP+ESUM< I )*10 . 

42 CONTINUE 

IF<LFLG<4> .EQ. 1) EXP=-EXP 

43 DO 50 I=1,K 
J=L+1-I 

IF<I .GT.L> J=I 

IF<LFLG<1 > .EQ. 1 > GO TO 45 

J=K+1-I 

INUM=INUM+TSUf-K I )*SN< J) 

GO TO 50 
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AN6 


67 

70 

71 

72 

73 

74 

75 

76 

77 
IOC 
101 
102 

103 

104 

105 

106 
107 
110 
111 
112 

113 

114 

115 

116 
117 

-.120 

121 

122 

*123 

124 

125 

126 


45 RNUti = RNuM + TSUM< I >>»«SN<:j)flO . **EXP 
50 CONTINUE 
C Ai .6 

C TEST FOR NEGATIVE INPUT NUMBER 
C 

IF<LFLGC1> ,EQ. 1) GO TO 53 
IF(LFLGC2) . EQ . n INUM=-INUM 
IINUM=INUM 
GO TO 55 

53 rF<LFLG<2) .EQ-. 1> RNUM=-RNUM 
RRHUM=RNUM 

55 IF<LFLG<8) ,EQ.O> RETURN 
LFLG<3>=1 

56 RETURN 

60 IFCLFLGCS) .NE. 1 ) GO TO 61 
LFLG<3)=1 

RETURN 

61 IFCICNT.GE. 10> RETURN 

IF<<LFLG<7> .EQ.O) ,XOR. <LFLGC2> .EQ.O>> GO TO 62 
RETURN 

62 ICNT=ICNT+1 
KCNT=ICNT 
J=1 

IF<ICNT.LE.5) GO TO 63 

J=2 

KCNT=ICNT-5 

63 NDX=6t'CKCNT-l > 

NME(J)=<NME< J> . AND. MASK<KCNT)> . OR . < , NOT . MASK <KCNT> .AND. 
* (LETTER .R.NDX>> 

LFLG<7>=0 

LFLG<2>=0 

RETURN 
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1 SUBROUTINE ANRST 

2 C A2.0 ^ . 

3 C RESET GRID INPUT VARIABLES 

4 C 


5 


K*L=M=II=ICNT=I inum=o 

6 


RRNUM=0 . 

7 


NME<1)=NME<2)=' 

lo- 


DO 60 1=1/10 

ll 


LFLG<n=0 

12 

60 

CONTINUE 

13 


DO 70 1=1/20 

14 


ESUM<I)=0 . 

15 


TSUM<I>=0 . 

16 

70 

CONTINUE 

17 


RETURN 

20 


END 



28 


VER 0 REV 


16 JUN 76 


PAGE 1 


SUBROUTINE AN? 

Nl.O ■ * 

CONSTRUCTION OF ALPHA/NUMERIC INPUT FOR UABCOM 

COMHON/'PARA/ LFLG< 10 ) i MFLACi<4) , SUM. LETTER. NO, NFLC / IN0,NFLAGa2) 
* / JJJ,KKK,LLL/MMM.NNN.NAME<2.60),NME<2), IINUM.RRNUK 
COMMON/LTR/ I FLAG < 40 > 

DIMENSION TSUM<20) . SNC20 ) ,ESUM<20) . MASK<5> 

LOGICAL IFLAG 

DATA TSUM/20*0 ./'.SUM/'O ./',LFLG<6)/0*',LFLGC3)^1/ 

DATA SN/20*0 ./'.ESUM/20+0 JCNT/O/ 

DATA MASK/OB. 7700000000B,7777000000B,7777770000B,77777777008^ 
JCNT=JCNT+1 

IFCJCNT.EQ. 1> GO TO 100 
IF<LFLG<3) .EQ.O) GO TO 89 
100 CALL ANRST 


VER 0 REV 16 JUN 76 


1 C 

2 C 

3 C 

4 

5 C 

6 C 

7 C 
10 

11 c 

12 C 

13 C 

14 

15 C 

16 C 

17 C 
20 

21 C 

22 C 

23 C 

24 


25 

20 

lFaFLG<5> .EQ 

. 1 > 

GO 

TO 

35 

26 


K=K+1 





27 


IF<LFLG< 1 > ,NE 

1) 

GO 

TO 

25 

30 

22 

ICNT=ICNT+1 





31 


II=1-ICNT 





32 


L=K+II 





33 


IF<ICNT.EQ. n 

GO 

TO 

24 


34 


GO TO 30 





- 35 

24 

LFLG<7)-0 





36 


GO TO 999 





37 

25 

II=K-1 





40 

30 

TSUM<K)=SUM 





41 


SN 00 = 10 .**11 





42 


LFLG<7)=0 





43 


GO TO 999 





44 

35 

^f=H+l 





45 


ESUM<M)=SUM 





46 


LFLG<7>=0 





47 


GO TO 999 





50 

999 

EXP=INUM=0 





51 


RNUM=0 . 





52 


IF<LFLG<5:) . NE . 

1 > 

GO 

TO 

43 

53 


DO 42 1=1. M 





54 


J=M-I 





55 


EXP=EXP+ESUM< I 


0 . * 

*<J 


56 

42 

CONTINUE 





57 


IF<LFLG(:4> .EQ . 

1> 

EXP 

=-EXP 

60 

43 

DO 50 1=1. K 





61 


J=L+1-I 





62 


tF< I .GT .L> J=I 





63 


IF<LFLG<1> .EQ. 

1 > 

GO 

TO 

45 

64 


J=K+l-I 





65 


INUM=INUM+TSUM 

<I) 

*SN 

«C J) 


66 • 

- 

GO TO 50 






N1 . 1 

ONE DATA ENTRY PER LIGHT PEN PICK 

89 IF<LrLG<6) .EQ. 1) GO TO 89 
N1 .2 

WHEN CLR PICKED BY LIGHT PEN 

1F<LFLG(9> .EQ. 1) GO TO 100 
Nl .3 

CONSTRUCTION OF ALPHA WORD 

irUFLAGOD) GO TO 60 
Nl ,4 

WHEN ENT IS PICKED BY LIGHT PEN 

IFCLFLG<8> ,EQ. n GO TO 999 
Nl .5 

DETERMINE IF NUMERICAL INPUT IS REAL OR INTEGER 
IF<LFLG<75) 22.56,20 
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67 

70 

71 

72 

73 

74 

75 

76 

77 
100 
101 
102 

103 

104 

105 

106 
107 
110 
111 
112 

113 

114 

115 

116 
117 

.-120 

121 

122 

*123 

124 

125 

126 


45 RNIJM-RNIJM+TSUMCI>'^SN<J>'»'10 .hc^EXP 
50 CONTINUE 
C N1 .6 

C TEST FOR NEGATIVE INPUT NUMBER 
C 

IF<LFLG<n .EQ. 1 > GO TO 53 
IF<LFLG<2) .EQ. n INUM=-INUM 
IINUM=INUM 
GO TO 55 

53 IFCLFLG<2) .EQ. 1 ) RNUM=-RNUM 
RRNUM=RNUM 

55 IF(LFLG(S) ,EQ.O) RETURN 
LFLG<3)=1 

56 RETURN 

60 IF<LFLG<8>.NE. 1> GO TO 61 
LFLG<3)=1 

RETURN 

61 IFCICNT . GE . 10 ) RETURN 

IF< <LFLG<7) .EQ.O ) .XOR. <LFLG<2> . EQ.0)> GO TO 62 
RETURN 

62 ICNT=ICNT+l 
KCNT=ICNT 
J=1 

IF<ICNT.LE.5> GO TO 63 
J=2 

KCNT=ICNT-5 

63 NDX = 6+-'CKCNT-1 ) 

NME< J:> = (NME< J> . AND.MASK(KCNT)) . OR . < . NOT . MASK<KCNT> . AND . 
* (LETTER. R.NDX>) 

LFLG<7>=0 

LFLG<2)=0 

RETURN 
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I SUBROUTINE fiHRST 


2 

C 

N1 .7 

3 

C 

RESET GRID INPUT VARIABLES 

4 

5, 

c 

K*L*M= 1 1 = I CNT= 1 1 NUM»0 

6 


RRNUM=0 . 

7 . 


NME< 1 )=NME<2)=' ' 

10 


DO 60 1=1. 10 

11 


LFLG<I>=0 

12 


60 CONTINUE 

13 


DO 70 1=1,20 

14 


ESUH<O = 0. 

15 


TSUM<I)=0 . 

16 


70 CONTINUE 

17 


RETURN 

20 


END 
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PACE 4 


AN7 VER 0 REV 


1 

2 

3 

4 

5 

6 
7 
10 
11 
12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

55 

56 

57 
60 
61 
62 

63 

64 

65 

66 


IMAGE 6RIDU 
C R 1 . 0 

C DISPLAY ALPHA/NUMERIC INPUT GRID 
C 

COMMON/PARA/' LFLG< 10 > > MFLAG<:4> , SUM. LETTER. NO. NFLG. INO. NFLAGC 12> 

* . JJJ.KKK.LLL.MMM.NHH.NAME<2.60>.NMEC2>. I INUM.RRNUM 
COMMON/GRD/' BBOX<50 TEST <200 J . TEST3CS0 >. OPDPL ClOO ). BOX<20 > 
INTEGER ALPHA<27).iPtlJD4 

DATA ALPHA-'IHA. IHC. IHM. IHS. 1H8. IHH. IHN. IHT. IHC. IHI. IHO.IHU. 

* IHD. IHJ. IHP. IHV. IHE. IHK. IHQ. IHU. IHY. IHF. IHL. IHR. IHX. IHZ. IH^/' 
DATA LFLG<10).'0/ 

LDXCl .F) 

LDY<-1 .F) 

LSCL<.4F> 

TABLE2D<B80X> 

CHAR<TEST3) 

IF<LFLG<6> .EQ.O) GO TO 200 
IF<<^LPNSHJ.AND.200B) .NE.O) GO TO 200 
LFLG<6)=0 
200 CLEARPEN<1) 

♦PWD4=0 
PENOH< 1 > 

MEHU< ‘‘QBaXQ0 61SYQ0 70 A9 . S A5X3 0 66aYaO 70 G8 , OAQ 

* 8X90739Y8C!?0Ma . 9A3X90 773Y90 70 S3 . aA3 

* 8X80 613Y3 0 74B9 , 9A9X30 66QYaO 74H3 . 3 A3 

* SX80739Y3074N3. 3A3X30773Y3074T3.3A3 

* SX30 6 1 SYS 100 C9 . 3A3X30 663Y3 1 0 0 1 3 . 3A3 

* 8X30733Y310003 , 3A3X30 773 Y3 1 0 0 U3 . 3 A3 

* 8X80613Y3106D3 . 3A3X30 663 Y3 1 0 6 J3 . 3 A3 

* 8X30 733 Y3106P3 . 3 A3X30 773Y3 1 0 6V3 . 3 A3 

* SX30613Y31 12E3 . 3A3X30663Y31 12K3 , 3 A3 

* 3X30 733 YSl 12Q3 . 3A3X30 773 Y3 1 1 21U3 , 3 A3 

* 3X31043Y31 12Y3 . 3A3X30 6 1 3 Y3 1 20 F3 , 3 AS 

* 8X80663Y3120L3.3A3X30733Y3120R3.3A3 

* 3X30 778Y8 120X3 . 3A3X3 1 0 43 Y3 1 20 23 . 3A3 

* 3X81 1 l8Y3l20/'3 . 3A8 

* 8X31043Y3070 13 . 3A3X31 1 13Y307023 . 3A3 

* 8X31153Y307033. 3 A3X3 1 0 43 Y30 7443 . 8 A3 

* 8X8 lIl3Ya07453.3A3X3l 153Y307463 , 3 A 3 

* 8X81043Y910073.3A3X31113Y310083.3A3 

* 3X81 153Y310093 . 3A8X31048Y310603. 3A3 

* 8X81 1 13Y3106*3 . 3A3S3X31 153Y3106CLR3 . 3A8 
>*« 8E3X81 1 13Y31 12-3 . 8A3S3X8 1 1 68 Y3 1 I 2EX3 . 3A8 

* 8S3X31 153Y3120ENT8 . 32“ ) 

PENOFFU) 

IF<PENHIT<1)> LFLG<10)=1 

GO TO < 16. 16. 16. 16. 16. 16. 16. 16. 16. 16. 16. 16, 16.- 16. 16. 16. 

16. 16. 16. 16. 16. 16. 16. 16. 16. 16. 16. 

* 1. 2. 3. 4. 5. 6. 7. 3. 9. 10. 11. 12. 13. 14. 15 > 

RETURN 

1 IF<LFLG<6) .EQ. 1> RETURN 
SUM= 1 . 

LETTER=1H1 

LFLG<6)=1 


<^N7 
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67 

LPLG<7)=1 


70 

♦ Pli)D4=0 


71 

RETURN 


72 

2 IF<LFLG<6) .EQ. li> 

RETURN 

73 

SUM =2 . 


74 

LETTER=1H2 


75 

LFLG(6>=1 


76 

LFLG<7>=1 


77 

♦ PUJD4 = 0 


100 

RETURN 


101 

3 IF<LFLG<6> .EQ. 1> 

RETURN 

102 

SUM=3. 


103 

LE7TER=1H3 


104 

LFLG<6)=1 


10 5 

LFLG<7>-1 


106 

♦ P(JiD4=0 


107 

RETURN 


110 

4 IF<LFLG<6) . EQ. 1 > 

RETURN 

111 

SUM=4 . 


112 

LETTER=1H4 


113 

LFLG<6)=1 


114 

LFLGC7>=1 


1 15 

♦ PUID4 = 0 


116 

RETURN 


117 

5 IF<LFLGc:6> .EQ. 1) 

RETURN 

120 

SUM=5. 


121 

LE7TER=1H5 


122 

LFLG<6)=1 


' 123 

LFLG<7)=1 


124 

♦ PW[>4=0 


125 

RETURN 


126 

6 IF<LFLG<65 . EQ. 1 ) 

RETURN 

127 

SUM=6. 


130 

LETTER=1H6 


131 

LFLG<6)=1 


132 

LFLG<7)=1 


133 

♦PUJD4=0 


134 

RETURN 


135 

7 IFCLFLG<6) . EQ . 1 ) 

RETURN 

136 

SUM=7. 


137 

LET7ER=1H7 


140 

LFLG<6)=1 


141 

LFLG<7>=1 


142 

iPUJD4=0 


143 

RETURN 


144 

8 IF<LFLG<6) . EQ . 1 > 

RETURN 

145 

SUM=S. 


146 

LETTER=1H8 


147 

LFLG<6)=1 


150 

LFLG<7)=1 


151 

♦PUD4=0 


152 

RETURN 


153 

9 IF<LFLG<6> . EQ. 1 ) 

RETURN 

154 

SUM=9. 
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155 


LETTER=1H9 

156 


LFLG<6)=1 

157 


LFLG<7)=1 

160 


♦PUD4=0 

161 


RETURN 

162 

10 

IFCLFLGC6) 

163 


SUM=0 . 

164 


LETTER=1H0 

165 


LFLG<6)=1 

166 


LFLG<7)=1 

167 


iPWD4=0 

170 


RETURN 

171 

11 

IF<LFLG<6) , 

172 


LFLG<1 )=l 

173 


LFLG<:6) = 1 

174 


LFLG<7)=-1 

175 


LETTER^IH-^ 

176 


♦PWD4=0 

177 


RETURN 

200 

12 

IF<LFLG<6) . 

201 


LFLG<6)=1 

202 


LFLG<9)=1 

203 


tPUJ04 = 0 

20 4 


RETURN 

205 

13 

IF<LFLG<6) , 

206 


IF<LFLG<5) . 

20 7 


LFLG<6)=1 

210 


LFLG<4)=1 

21 1 


LFLG<7)=0 

212 


iPW04=0 

213 


RETU.RN 

214 

18 

LFLG<2)=1 

215 


LFLG<7)=0 

216 


LFI.G<6) = 1 

217 


LETTER=1H- 

220 


♦PWD4=0 

221 


RETURN 

222 

14 

if<:lflg<6) . 

223 


LFLG<6)=1 

224 


LFLG<5)=1 

225 


LFLG<7)=0 

226 


iPljJD4=0 

227 


RETURN 

230 

15 

IFCLFLGCE) .1 

231 


LFLG<6)=1 

232 


LFLG<8)=1 

233 


$PWD4=0 

234 

■ 

RETURN 

235 

16 

IF<LFLG<6) .i 


16 JUN ?6 


236 

237 

240 

241 

242 


EQ.l) RETURN 


EQ.n RETURN 


EQ,1) RETURN 


EQ.l) 
NE. 1 ) 


RETURN 
GO TO 18 


EQ.l) RETURN 


EQ.l) RETURN 


EQ.l) RETURN 


LFLG<6)=1 

LFLG<7)=1 

LETTER=8LPHA<tPWD4) 

♦PW04=0 

RETURN 


PAGE 4 
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END 
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PAGE i 


1 

2 C 

3 C 

4 C 

5 

6 
7 

10 C 

11 c 

12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 10 

26 20 
27 

30 

31 

32 


SUBROUTINE BLDNAM<NAM> 

J1 .0 

FORMATS FILE INPUT NAME FOR FILE ACCESSING SUBROUTINES 

IMPLICIT INTEGER <A-2) 

DIMENSION NAM<3) 

DATA MASK/778/ 

IF<NAM<2) . EQ . IH > NAM<2>=NAM< 1 ) 

IF<NAM<1>.EQ.NAM<:25) RETURN 
NAM<3)=NAM<2) 

NAM<2)=NAM< 1 > 

DO 10 I=9.0.-l 

UID*I/5+2 

KNT=<6*I >^/'30 

ISHf T=24~KNT 

TM=MASK. L . ISHFT 

POSC=NAM<WD> .AND.TM 

CHAR=CPOSC.L.KNT) .OR. IHC 

IF<<CHAR.NE. IHC) .AND. (CHAR.NE, IH >> GOTO 20 
NAM<WD)=NAM<WO) , AND . <-TM) 

CONTINUE 

CONTINUE 

NAM<l) = ONAM+l) .AND.77777B 
NAM< 1 )=NAM< 1 > .OR.7777500000B 
RETURN 
END 
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VER 0 
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1 SUBROUT I ME BLDNME<NAME> 

2 C Bl.O 

3 C FORMATS FILE INPUT NAME FOR FILE ACCESSING SUBROUTINES 

A C 

5 DIMENSION L< 10 > . NAME<3) 

6 NINE=9 

7 READCNINE. 100> L 


10 

100 

FORMAT<10R1) 

11 


DO 3 M=1.3 

12 

3 

NAMECM)=0 

13 


J=2 

14 


NL=0 

15 


NBLNK=0 

16 


DO 50 1=1.10 

17 


IF<L< I > . EQ. 40B) GO TO 15 

20 


IF<NBLNK.NE.O) GO TO 20 

21 

5 

IDUM=L< I > 

22 


CALL IT 

23 


GO TO 50 

24 

15 

NBLHK=NBLNK+1 

25 


GO TO 50 

26 

20 

IDUM=40B 

27 


DO 21 L=1.NBLNK 

30 

21 

CALL IT 

31 


NBLNK=0 

32 


GO TO 5 

33 

50 

CONTINUE 

34 


IF<NBLNK,EQ.O) GO TO 60 

' 35 


IDUM=0 

36 


DO 70 L=1.NBLNK 

37 

70 

CALL IT 

40 


_ HBLHK=0 

41 

60 

IF<NAMEf3) . NE . 0 > GO TO 80 

42 


NAME< n=NAME<2) 

43 


RETURN 

44 

80 

L0C:=aNAMEi:2) 

45 


NAME<n = < .H.77775B) .0. LOC 

46 


RETURN 

47 

C 


50 


SUBROUTINE IT 

51 


NAME< J>= ,K .NAME< J> 

52 


NAME< J>=NAME< J) ,0. IDUM 

53 


NL=NL+1 

54 


IF<NL.GE,5) J=3 

55 


RETURN 

56 


END 
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3 

4 

5 

6 
7 
10 
11 
12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 
'35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

55 
' 56 

57 

60 

61 

62 

63 

64 

65 

66 


SUBROUTINE CGMOVE ^ 

Cl.O 

BLACK BOX TRACKER CG TRANSLATOR 

C0MM0N/'M0VIT/NAME<2) >XCG. YCG.ZCGi NREAD 

COMMON/'RSTORE/CAMB. ALEN/ RAD2<30 > , RADY <30 ) . RADXOO ) , CAMX<30 ) , ’ 

CAM2<30 > . FACTOR. I BLOCK . F I LLER < 56 ) 

LOGICAL FLAG. ADTF1-. ADTF2 . TOP . NREAD 
LOGICAL FPEN. FPENP. NDSTP 

FRACTION RADX. RAD2. SX. SY. RADTX. RADTY . POSTX , POSTY. P0ST2 
FRACTION RX. RY. RADY. BIAS. FDUMOO ) 

DIMENSION NAMDUM<3) . IDUMOO ) . NOTES <50 ) . NOTES 1 < 50 ) . N0TES2 < 50 > 
DIMENSION TVIEW< 100 ) . SV I EUl < 1 0 0 ) . DUMP 1 T < 2 1 0 ) 

EQU I VALENCE< DUMP I T . CAMB > 

EQUI VALENCE<NPTS. RAD2< 1 ) ) . <NPTSC. CAM2< 1 ) ) . <NPTSY. RADY< I >) 


♦NDSTP=. TRUE. 

Cl . 1 

DINGER DATA 

DATA2D-<DINGER> 
2SET< 1 ,F> 

M0VE<1 ,F. 1 .F) 
DRA0J<1 .F.-l ,F> 
DRAU<-1 .F. -1 ,F) 
DRAtiK-l ,F. 1 .F> 
DRAU< 1 ,F. I . F> 
ENDLIST 
ENDDATA 


IF<NDSTP> GO TO 1 
IF<NREAD> GO TO 2000 
Cl .2 

INITIALIZE CG POSITION 

POSTX=-RTDF<XCG/FACTOR+ .8) 

POSTY=RTOF<YCG.'FACTOR) 

P0ST2=RT0F<ZCG/FACT0R)- 
Cl .3 

INITIALIZE LABELS 

LABELCNOTES) 

URITE< 16. 40 ) NAME 

FORMAT <"863X30 30 3Y30 14COMPONENTSXa074XCGaX8l 15YCG9X81 352CG8XS0 30 SY 
*S020".2A5. "3. ") 

ENDLIST 
LABEL < NOTES!) 

WRITE < 16. 41 ) XCG. 2CG 

rORMAT<"3BsXS065aYs020".F7. 1. "SXS125".F7. 1. "8. ") 

ENDLIST 


39 


CGMOVC 


VCR 0 REV 
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pr^GE 1 


67 


LABELCNOTES2) 

70 


UIRITE<16,42)YCG 

^71 

42 

^ FORMAT<"QB8Xal05SY3020’', F7. 1/ "8. 

72 


endXT^t 

73 

86 

DISPLAY TWERP 

74 

1000 IFCDCPRESSC 1 > > GO TO 6000 

75 

C 

Cl .4 

76 

c 

- WHERE IS DINGER? 

77 

c 


100 


RX=4RADTX 

101 


RY=4RADTY 

102 


TOP=. FALSE. 

103 


IFCRADTY.GT.O .F) TOP=,TRUE. 

104 


XPOS=FTOR<RX)+ . 8 

105 


BIAS= . 6F 

106 


IFCTOP) BIAS=-.4F 

107 


YPOS=FTOR<RY+BIAS) 

110 

c 

Cl . 5 

111 

c 

DINGER POSITION IN ENGINEERING UNI 

112 

c 


113 


ENGX=-XPOS*FACTOR 

114 


ENG2=YP0S*FACT0R 

115 


DIMENSION POSNaO) 

116 

c 

Cl .6 

117 

c 

DINGER POSITION ON SCREEN 

120 

c 


121 


LABEL<POSN> 

122 


IF<TOP> GO TO 36 

*123 


WRITEa6,89) ENGX.ENGZ 

124 

89 

FORMAT< ''8S0BaXaO673Y9l5OX=" , F6 . 0, 

125 


GO TO 37 

126 

36 

WRITEC16>91) ENGX.ENGZ 

127 

91 

FORMATC ■•8SaBaXS0 67aYal50X=". F6. 0. 

130 

37 

CONTINUE 

131 


ENDLIST 

132 


IF< .N.FLAG> GO TO 1000 

133 


NDSTP= , TRUE. 

134 

200 

IF<NDSTP) GO TO 200 

135 

C 

Cl .7 

136 

C 

MOVE CG TO NEW POSITION 

137 

C 


140 


IF<TOP) GO TO 210 

141 


POSTY=SY+. 6F 

142 


YCG=FTOR<POSTY)+FACTOR 

143 


LABELCN0TES2) 

144 


WRITE< 16, 42)YCG 

145 


ENDLIST 

146 


FLAG=. FALSE. 

147 


GO TO 86 

150 

210 

POSTX=SX 

151 


P0ST2=SY-. 4F 

152 


XCG=-<FT0R<P0STX)+ . 8)*FACT0R 

153 


2CG=FT0RCP0ST2)*FACT0R 

154 

- 

LABELCNOTESl ) 


=“.F6.0> 

=".F6.0> 
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PR6E 1 


155 


URITE<16, 41 )XCG.2CG 

156 


ENDLIST 

157 


FLAG=. FALSE. 

160 


GO TO 86 

161 

C 


162 

6000 

CONTINUE 

163 


RETURN 
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PRGC 2 


1 

2000 

CALL iPTCH^ 

2 

C 

Cl .8 

3 

C 

GET FUSELAGE SHAPE FROM DISK 

4 

C 


5 


NINE=9 

6 


UJRITE<NINE,2001 ) 

7 

2001 

FORMAT<’ ENTER PACK/'VOL NO OF INPUT.. 

10 


READ-CNINE^O) LPV 

11 


LLPV=IOCTCLPV) 

12 

2015 

tDlRITE<NINE, 2002) 

13 

20 0 2 

FORMAT<‘ ENTER NAME OF INPUT FILE..’) 

14 


CALL BLDNME<NAMDUM) 

15 


$NEUIID = 5000000000B 

16 


♦DFILETYPE=17B 

17 


IFILAD=1 

20 


CALL RDSKOKCLLPV, NAMDUM. 208, IFILAD, IFLG) 

21 


IFdFLG.NE.O) GO TO 2010 

22 


HIRITE< 10 , 20 03)NAMDUM(2) , NAMDUM < 3 ) , LLPV 

23 

2003 

FORMAT<’ CAN NOT LOCATE ',2A5,' ON PACK/VOL 

24 


GO TO 2010 

25 

2010 

CALL RNDRD<DUMPIT, 0,208, IFILAD, lOFF) 

26 


CALL DSKER 

27 


NEWID=7000000000B 

30 


DFILETYPE=7B 

31 


CALL SKETCH 

32 


NREAD=, FALSE. 

33 


GO TO 72 

34 

C 



,04) 
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1 

C 


2 

C 1 ! 

I ! > t ! : ! : t ! 1 : i : t ! t ! ; t t i t t ; t ] 1 t 1 : i ! i 1 : I { 1 : t > 1 1 t : i 

3 

C 


4 


iUt^ftOUlINE SKETCH 

5 

C 

Cl .9 

6 

C 

CREATE FUSELAGE OUTLINE 

7 

C 


10 


ISS=NPTS+1 

11 


DO 2 IS=2> ISS 

12 


FDUMaS)=-RAD2aS) 

13 

2 

IF<CAMB> FDUMaS)=FDUMaS>+CAMZ<IS>+CAMZaS) 

14 


LIST2D<SVIE!lO 

15 


ZSET<0 . F) 

16 


DO 3 IS=2, ISS 

17 


IF<IS.GT.2> GO TO 3 

20 


MOVE<RADX< IS). RADZ< IS> ) 

21 

3 

DRAUJ<RADX<: IS) . RADZCIS)) 

22 


DO 4 IS=2. ISS 

23 


IF<IS,GT.2) GO TO 4 

24 


MOVECRADX< IS) .FDUMCIS)) 

25 

4 

DRAUKRADX< IS) .FDUM<IS>) 

26 


EHDLIST 

27 


LIST2D<TVIEUJ) 

30 


ZSET<O.F) 

31 


DO 5 IS=2. ISS 

32 


IF<IS,GT.2) GO TO 5 

33 


MOVECRADXCIS) . RADY < IS) ) 

34 

5 

DRAUKRADX<IS).RADYaS)) 

35 


DO 6 IS=2. ISS 

36 


IFCIS.GT.Z) GO TO 6 

37 


MOVE<RADX< IS). -RADYC IS)) 

40 

■9 

DRHliKRADX<IS). -RADY< IS)) 

41 


■ ENDLIST 

42 


RETURN 
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1 

C 


2 

Cl ; 


3 

C 


4 


IMAGE TWERP 

5 

c 

C2.0 

6 

c 

DISPLAY FUSELAGE 

7 

c 


10 


READ<3> 

11 

c 

C2. 1 

12 

c 

CHECK DDT PEN FOR TIP SWITCH 

13 

c 


14 


FPEN=. FALSE. 

15 


IF<,N.4ADTFi .A. $ADTF2) FPEN=. 

16 


IF<FPEN.A. .N.FPENP) FLAG=,TRUE. 

17 


IFC .N. FPEN) GO TO 5 

20 

c 

C2.2 

21 

c 

SAVE DDT PEN POSITION 

22 

c 


23 


SX=tRADTX 

24 


SY=iRADTY 

25 

5 

FPENP=FPEN 

26 

c 

C2.3 

27 

c 

GRID DATA 

30 

c 


31 


DATA2D(GRD> 

32 


ZSET<1.F> 

33 


LINE<-.8F.O .F, .82F.0 .F) 

34 


LINEC-.8F.0 ,F,-.8F, .41F) 

‘ 35 


LINE<-. 8F. 0 ,F, ,82F. 0 . F> 

36 


LINE<-.8F.O .F,-.8F, .41F) 

37 


ENDLIST 

40 


ENDDATA 

41 


' LDYC . 4F> 

42 

c 

C2. 4 

43 

c 

SHOW SIDEVIEW & GRID 

44 

c 


45 


TABLE2D<SVIEW> 

46 


TABLE2D<GRD) 

47 


P0SCHAR<-,82F. .42F/0 .F. "9S3B-Z" 

50 


POSCHARC .83F>-.01F,0 .F> ''sSoB-X" 

51 


P0SCHAR<-.2F> .45F.0 .F, "9BSIDEVI 

52 

c 

C2.5 

53 

c 

POSITION 3< DISPLAY SIDEVIEW CG 

54 

c 


55 


DY<POSTZ> 

56 


DX<POSTX) 

57 


LSCL< .IF) 

60 


TABLE2D<DINGER> 

61 

c 

C2.6 

62 

c 

SHOW TOPVIEW & GRID 

63 

c 


64 


LSCLXl.F) 

65 


LDY<- . 6F) 

66 

- 

LDX<0 .F) 
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PR6C 4 


67 

70 

71 

72 

73 

74 C 

75 C 

76 . C 

77 
100 
101 
102 

103 

104 

105 C 

106 C 

107 C 

no 

111 

112 

113 

114 

115 C 

116 C 

117 C 
■ , 120 

’ 121 
122 
■ 4 23 

124 

125 


TAeLE2D<TVIEW> 

TABLE2D<GRD> 

P0SCHAR<-.8F, .42F,0 ,F, "SSSBY") 
POSCHAR< ,82F.-.0 0 5F,0 .F, "seaS-X'^) 
POSCHARC - .2F> ,45F.O .F, " 8BPLANV I EU) " ) 
C2.7 

POSITION 8, DISPLAY TOPVIEW CG 

DXCPOSTX) 

DY<POSTY) 

LSCL< . IF) 

TABLE2D<DINGER) 

LSCLCl .F) 

CENTER 

C2.8 

WRITE NOTES ON SCREEN 

CHARCPOSM) 

CHARCN0TE31 ) 

CHAR<N0TES2) 

CHARCNOTES) 

IFC$ADTF1) RETURN 
C2.9 

DISPLAY DINGER 

LDXCRX) 

LDYXRY) 

LSCLC .OIF) 

TA8LE2D<DINGER) 

RETURN 

END 
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CROSS 


VER 1 
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1 SUBROUTINE CROSS < Ai B . C, CM> 

2 C 

3 C UNITIZED MATRIX CROSS PRODUCT C « A X B 

4 C 

DIMENSION A<3>,B<3).C<3) 

C< 1 )=A<2)*8<3)-A<3)'*-B<2) 

C<2)=A<3)*BC1 )“Aa )*B<3) 

10 C<3) = A<1>*B<2)“A<2)>*-Ba) 

11 CM=SQRT<CC1 >*C<1 >+CC2>>*‘C<2)+C<3>fC<3)> 

12 RETURN 

13 END 
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1 SUBROUTINE DEC FMT <VALUE. NFIELD. FMT> 

2 C Dl.O 

3 C BUILD FORMAT FOR DECIMAL NUMBERS 

4 C 

5 INTEGER FMT. CVT<2) . POINT 

6 DATA <CVT<I>. 1 = 1.2) .'IHF.IHE^ 

7 . DATA POINT/IH./ 

10 C * ♦ DETERMINE THE EXPONENT <NPUR 

11 NPUR=0 

12 ABVAL = A8S<< 1 . +10 . **<3-NFIELD) )>*«VALUE) 

13 IF <ABVAL.LT. 1 . ) GO TO 550 

14 DO 530 1=1.99 

15 ABVAL=ABVAL/'lO . 

16 IF <ABVAL.LT. 1 . > GO TO 540 

17 530 CONTINUE 

20 540 CONTINUE 

21 NPUR=I 

22 GO TO 600 

23 550 CONTINUE 

24 DO 570 1=1.99 

25 ABVAL=ABVAL+10 . 

26 IF CABVAL.GT. . 1) GO TO 580 

27 570 CONTINUE 

30 58Q CONTINUE 

31 NPUJR = -I 

.'-32 C ♦ * SELECT THE FORMAT 

33 600 CONTINUE 

-34 NSP=1 

35, IF <VALUE.LT.O) NSP=2 

36 IF <NPUR,GE.O) GO TO 620 

37 IF <NPlJiR.LT.-4) GO TO 640 

40 • NDEC=NFIELD-NSP 

41 IFMT=1 

42 SO TO 690 

43 620 CONTINUE 

44 NDEC=NFIELD-NPWR-NSP 

45 IF <NDEC.LE.O) GO TO 640 

46 IFMT=1 

47 CO TO 690 

50 640 CONTINUE 

51 NDEC=NFIELD-NSP-5 

52 IFMT=2 

53 690 CONTINUE 

54 IF CNDEC.LT.O) NDEC=0 

• 55 C ★ ♦ BUILD THE FORMAT 

56 ENCODE <5. 1000. FMT) CVTC IFMT) . NFIELD. POINT. NDEC 

57 1000 FORMAT < A1 . 12. A 1 . 1 1 ) 

60 RETURN 

61 END ' 
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1 

2 

3 

A 

5 

6 
7 
10 
11 
12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

55 

56 

57 
60 
61 
62 

63 

64 

65 

66 


image grid 
C G1 .0 

C DISPLAY ALPHA/'NUMERIC INPUT GRID 
C 

COMMON/'PARA-' LTLG< 10) / MFLAG<4) , SUM. LETTER / NO, NPLG / INO, NFLAG< 12> 
m , JJJ,KKK,LLL,MMM,NNN,NAME<2,60),NME<2), IINUM.RRHUM 
COMMOH/GRD,' BBOX < 50 ), TEST , TESTS 
INTEGER AL'PHA<27), 4PWD4 

DATA ALPHA/IHA, IHG, IHM, IHS, IHB, IHH, IHN, IHT, IHC, IHI, IHO, IHU, 

* IHD, IHJ, IHP, IHV, IHE, IHK, IHQ/ IHU, IHY, IHF, IHL/ IHR, IHX, IHZ, 1H>'/ 
DATA LFLG<;i0)/0^ 

LDXO .F>- 
LDY<-1 . F) 

LSCL< .4F) 

TABLE2D<BB0X) 

CHAR<TEST) 

IF<LFLG<6) .EQ.O) GO TO 200 
IF<<iLPNSW.AND.200B) .NE.O) GO TO 200 
LFLG<6)=0 
200 CLEARPEN<1> 
iPWD4=0 
PENOfKl) 

MENU< "3B5iXa0 6l6YQ0 70A9 , 3 A3X80 66sYg0 70 Gs .8 AS 

* 8X80733Ya070Ma.SAaXs0778Y8070S8.SA9 

* SXs06laYa074B8.aAaXs0668Ys074H8.SAa 

* 8X80 73aY80 74Na . SAaXaO 778Y80 74T8 . 8 AS 

* 8X80 618Y8100C8 . SASXaO 668 Y3 1 0 0 IS , 8AS 

* 8X80 73SY810008 . aAsXao 77a Ya I 0 0 ua . a AS 

* sXa061aYal06Da . 8AaX8066aY8l06Ja , aAa 
>♦« 8X80733Y8i06PS . 8A8Xa0778Yal06V8 .8A8 

* 8Xa06l8Y8l 12E© 8AaXa066SY8l 12K8 , 8A3 
. >*« 8Xa0738Yall2Qa.8AaXa0778Y8ll2Ula.8A8 

* SXal04aYall2Y8.8AaX306l3Y8l20Fa,8Aa 

* SXSO 66SY8120L8 . 8AaXS0 738Ya 1 20 R3 . 8A8 
>♦' 8Xa077aYai20X8.3Aaxai04aY8i202a.3Aa 

* sXal 1 lsYai20/'a , aA8 

* sXalo 4aYa0 70 l a . sAaXal l laYa0 70 23 . aAa 

* 8X81 15aY80703a , aA8XSl04aY307448 , SAa 

* 8X81 1 laY307458 . 8A8X81 15aY807468 . 8A8 

* 8Xai04aYai007a . aAaXai i isYalooSa , sAa 

* 8X8li5aY3i009a.3A3Xai04aYai060a.aAs 

* 8Xai 1 laYsios+st . aAasaxal i3sYai0 6CLRa . 8A8 

* 8E8X81 1 laval 12-a . 3AaSaX8l 168Y81 12EX8 . 3A8 

* 8S8X8ll5aYal20ENTa.a2"> 

PEHOFF< 1 ) 

IF<PCNHIT< 1 ) ) LFLGU0) = 1 

GOTO <16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16, 

16, 16, 16, 16, 16, 16, 16, 16, 16,16,16, 

* 1,2,3,4,5,6,7,3,9,10,11,12,13,14,15) 

RETURN 

1 IF<LFLG<6).EQ.l) RETURN 
SUM=1. 

LETTER=1H1 

LFLGC6)=1 
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67 

70 

71 

72 

73 

74 

75 

76 


LFLG<7>«1 

iPUJD4=0 

RETURN 

!F<LFL6<6) . EQ . I 

SUM»2. 

LETTER=1H2 

LFLG<6)*1 

LFLG<7)=1 


77 


♦PUD4=0 

100 


RETURN 

101 

3 

IF<LFLGC6> .EQ. 1 

102 


SUM=3 . 

103 


LETTER=1H3 

104 


LFLG<6>=1 

105 


LFLG<7>=1 

106 


4PDJD4=0 

107 


RETURN 

110 

4 

IF<LFLG<6> .EQ. 1 

111 


SUM=4 . 

112 


LETTER=1H4 

113 


LFLG<6)=1 

114 


LFLG<7>=1 

115 


♦PWD4=0 

116 


RETURN 

117 

5 

IF<LFLGC6) .EQ. 1 

• 120 


SUM=5 . 

121 


LETTER=1H5 

122 


LFLG<6>=1 

' 123 


LFLG<7>=! 

124 


♦PUJD4 = 0 

125 


RETURN 

126 

6 

IF<LFLGC6) .EQ. 1 

127 


SUM=6. 

130 


LETTER=1H6 

131 


LFLG<6)=1 

132 


LFLG<7)=1 

133 


tPU04=0 

134 


RETURN 

135 

7 

IF<LFLG<6> .EQ. 1! 

136 


SUM=7. 

137 


LETTER=1H7 

140 


LFLG<6)-V- 

141 


LFLG<7>=1 

142 


iPlUD,4=0 

143 


RETURN 

144 

8 

IF<LFLG<6> ,EQ. I' 

145 


SUM=8 . 

146 


LETTER=1H8 

147 


LFLG(:6) = 1 

150 


LFLG<;7) = 1 

151 


<-PWD4=sO 

152 


RETURN 

153 

e 

IF<LFLG<6> .EQ. i; 

154 

- 

SUM=9 . 


> RETURN 


) RETURN 


) RETURN 


RETURN 


RETURN 


RETURN 


RETURN 


RETURN 
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155 

156 

157 
160 
161 
162 

163 

164 

165 

166 
167 

170 

171 

172 

173 

174 

175 

176 

177 
200 
20 1 
202 

203 

204 

205 

206 
207 
210 
211 
212 

213 

214 

215 

216 
217 
220 
221 
222 

223 

224 

225 

226 
227 

230 

231 

232 

233 

234 

235 

236 

237 

240 

241 

242 


LETTER=IH9 
l.FL/;f6>-l 
LFLG<7>=1 
♦PIU94 = 0 
RETURN 

10 IF<LFLG<6) .EQ. 1) RETUR?? 
SUM=0 , 

LETTER=1H0 

LFLG<6>=1 

LFLG<7>=1 

tPUD4=0 

RETURN 

11 IFCLFLG<6) ,EQ. 1) RETURN 
LFLG<1>=1 

LFLG<6)=1 

LFLG<7>s-l 

LETTER=1H* 

♦PWD4=0 

RETURN 

12 IFCLFLGCS) ,EQ, 1> RETURN 
LFLG<6>=:1 

LFLG<9)=1 

♦PWD4=0 

RETURN 

13 IF<LFLG<6> . EQ, 1 ) RETURN 
IF<LFLG<5) .NE. 1> GO TO 18 
LFLG<6)=1 

LFLG<4>=1 
LFLG<7)=0 
tP DID 4=0 
RETURN 
18 LFLG<2)=1 
■LFLG<7)=0 
LFLG<6)=1 
LETTER=1H- 
$PWD4=0 
RETURN 

14 IF<LFLG<6> .EQ. n RETURN 
LFLG<6>=1 

LFLG<5)=1 
LFLG<7)=0 
♦ PHJD4 = 0 
RETURN 

15 IF<LFLG<6) .EQ. 1) RETURN 
LFLG<6)=1 

LFLG<8>=1 

iPUD4=0 

RETURN 

16 IF<LFLG<6) . EQ . 1 :■ RETURN 
LFLG<6)=1 

LFLG<7>=1 

LETTER=PLPHA<iPUD4> 
tPWD4=0 
. RETURN 
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END 
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1 

2 

3 

4 


10 

11 

12 

13 

14 

15 


FUNCTION lOCKNN) 

C U1 .0 

C COMPUTES OCTAL EQUIVALENCE FROM DECIMAL INPUT 
C 

ITMP=0 

ITMP = <<NN/'1000 )*512) + ITMP 

NRMD=NN/^1000 

NTMP=< <NRMD^100)*64) + ITMP 

IlHP = NRMD/'/'100 

JTMP^C <ITMP/10)>4«8)+NTMP 

IOCT=<ITMP//10)+JTMP 

RETURN 

END 
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1 

2 

3 

4 

5 

6 
7 
10 
1 1 
12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

55 

56 

57 
60 
61 
62 
63 


IMAGE OPTIONS 
C 01.0 

C DISPLAY PROGRAM OPTIONS 
C 

COMMON/PARA/ LFLG< 10 MFLAG<4) / SUM. LETTER. NO. NFLG> INO. NFLAGC 12) , 
«*• JJJ.KKK. LLL.MMM,NNN.NAME<2.60).NMEC2). IINUM.RRNUM 
COMMON/OPT/' NHIT.OPDPL 
INTEGER iPUJD4 
DATA NHIT.'O/' 

DATA2D (BOX) 

ZSETCO.F) 

MOVE< .4F. .4F) 

DRAIIK .4F.-1 .F) 

DRAUJ<-.4F.-1 .F) 

DRAIi;<“,4F. .4F) 

DRAUK ,4F. .4F) 

ENDLIST 
ENDDATA 
LDXO .F) 

LDY<1.F) 

TABLE2D (BOX) 

CHAR< "3BaXQO?l3Ya0 640PTIONS'' ) 

CHAR(OPDPL) 

IFCNHIT.EQ.O) GO TO 100 

IF(<iLPNSW. AND.200B) .NE.O) GO TO 100 

NHIT=0 

100 CLEARPEH(l) 

$piUD4=0 

PENCN(l) 

MENU( "aeaSaX3061aY9072PACK/'VOLQ , aA8 
aX&061aYa076SELECT FILEa.SAS 
aXa0 6la?al0 2LIST VOLa.aAa 

* SXa061SYai06BLACK BOXa.aAa 

* axaoeiavai 12GE0METRY3 . aAa 

* aXa06laYai 16DISPLAY3 . aAa 

* aXa06laYai22APPEND DiSPLAYa.aAa 

>•« aXa061aVal26DELETE COMPONENTa . aA9 
t' 8Xa06iaYai32STOREa.aAa 

* 8X30 61 aval 36 INPUT NEUJ FILEa.aAs 

* 8Xa061aY3l42BLACK BOX TRACKERa . 8A8 

* 8Xa06laYai46PRlNTa . az" ) 

PENOFF(l) 

IF < <LF_G< 10 ) .EC. 1 ) . OR. (MFLAGd ) ,EQ . 1) .OR . <MFLAG<2) .EQ . 1) )GOTO 20 0 

IF(TPliJD4 , EQ . 0 ) GO TO 200 

IF<NHIT . EQ . 1 ) 60 TO 200 

NFLAG(’$.PIJIID4) = 1 

MFLAG(3)=1 

iPUJD4=0 

200 CALLSIMG GRID 
RETURN 
END 
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1 F INCLUDE=FORMT 

2 C Fl.O 

3 C FILE LOCATOR BY FILE SEQ. NUMBER IN PACK^VOL 

4 C 

5 SUBROUTINE LOCFLE < IPV/ NAMER. IFILNO) 

6 COMMONxPARA^ LFLG< 10 ) >MFLAG<4) > SUM, LETTER. NO. NFLG/ INO,NFLAG<12) 

7 . * .JJJ.KKK.LLL.MMM.NNN,HAME<2.60>,NME<2>. IINUM.RRNUM 

10 COMMON/'STR/' NTOT, value 1 <60 ) , VALUE2<60 ) , VALUES <60 > . VALUE4<60 ). 

U VALUE5<60).VALUE6<60>,VALUE7<60>,VALUE8<60)iVALUE9<60>. 

12 ♦ VALUEIO <60 > , VALUEl 1 <60 ) , VALUE12<60 ) . VALUE13C60 ) 

13 REAL VALUE1/VALUE2. VALUES, VALUE4, VALUES. 

14 * VALUES, VALUE?, VALUES, VALUES, VALUEIO. 

15 * VALUEl 1 , VALUE12, VALUE13 

16 INTEGER NAME, NAMER<3) 

17 CALL BLDNAM< HAMER) 

20 CALL DSVOL< IPV, IVFLAG.NFREES) 

21 CALL DSKER 

22 C URITE< 10 , 100 ) I VFLAG . HFREES 

23 C 100 F0RMAT< IX. 5013) 

24 CALL DSFIL<IPV,HAMER, IFFLAG.NMFIL. IFILAD. IFILNO) 

25 CALL DSKER 

26 C UIRITE<10,100)IFFLAG.NMF1L. IFILAD. IFILNO 

27 RETURN 

30 END 
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1 

2 

3 

4 

5 

6 
7 
10 
11 
12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

55 

56 

57 
60 
61 
62 

63 

64 

65 

66 


F ;nclude=formt 

SUBROUTINE INPUT<IPV> 

C II .0 

C INPUT NEU DATA FILE 
C 

COMMON/STR/ NTOT. VALUE 1 <60 ) , VALUE2<60 ) . VALUES <60 ) / VALUE4<60 > . 

♦ VALUE5<60>,VALUE6<60>,VALUE7<60),VALUE8<60>,VALUE9<60>, 

* VALUEIO <60 ) , VALUEl 1 <60 ) , VALUE 12 <60 ) . VALUE 1 3 < 60 > 

DIMENSION IBUF<1000).DUMN<2> 

REAL VALUEl . VALUE2> VALUE3 , VALUE4 , VALUE5 > 

♦ VALUE6. VALUE?, VALUES/ VALUES, VALUEIO. 

♦ VALUEl 1 / VALUE12/ VALUE13 
INTEGER NAME<2/60 )/NAMER<3) 

DATA DUMN<1)/'5H z' / DUMN < 2 ) /5H DUMX/'O . /'/ DUMY^O . DUMZ/'O . x 

DATA NINE/'?/' 

DATA FMT2/'l 51 51515156,^ / FMT3/7777777777B/ 

DATA FMT4/'????700000B/' 

NTOT=0 
CALL PTCH9 
UIRITE<NINE/ 3) 

3 FORMAT </'•/'/ ' ENTER NAME OF NEW FILE...') 

READ<NINE/5) NAMER 
5 F0RMAT<2A5) 

CALL BLDNAM< NAMER) 

10 1)JR1TE<NINE/ 1776) 

i.?6 FORMAT</// ' ENTER NAME OF COMPONENT...') 

NTOT=NTOT+l 

REACXNINE / 100) NAME < 1 . NTOT ) , NAME < 2 , NTOT ) 

100 F0RMAT<2A5) 

IF<NAME< 1 / NTOT) , EQ . 4HD0NE) GO TO 25 
WRITE<NINE/ 1976) 

1976 FORMAT<// ' ENTER 13 VALUES.,.') 

READ<NINE/ 0 ) VALUEl <NTOT) , VALUE2 <NTOT ) , VALUE3 < NTOT ) 

* / VALUE4<NT0T) / VALUE5<NT0T) , VALUE6<NT0T) / VALUE7<NT0T) 

* / VALUES < NTOT) / VALUES < NTOT ) , VALUE 1 0 < NTOT) / VALUE 1 1 < NTOT) 

* / VALUEl 2<NT0T) / VALUE 1 3 < NTOT ) 

GO TO 10 

25 NTOT=NTOT-l 

ENC0DE<5/ 1000 / IBUF<D) NT0T/FMT2 
-LOO F0RMAT<I3/A2) 

J=2 

DO 5000 1=1 /NTOT 

ENC0DE<75/ 2000 / IBUF<J)) NAME<1/ I)/NAME<2/ I ) > VALUEl < I ) / VALUE2< I ) / 

•i: VALUES < I ) / VALUE4< I ) , VALUE5< I ) / VALUES < I ) /FMT2 
2000 F0RMAT<2A5/ 6F10 . 1 / A5) 

ENCODE <75/ 30 00 / I EUF < J+ 1 5 ) ) VALUE? < I ) . VALUES < I )/ VALUES < I ) / 

♦ VALUElOa)/ VALUE1KI)/VALUE12<I)/VALUE13<I)/FMT2 
3000 FORMAT<7F10 . 1 / AS) 

ENC0DE<45/ 7000 / IBUF<J + 30)) DUMN< 1 ) / DUMN <2 ) / DUMX / DUMY / DUM2/ FMT2 
7000 F0RMAT<2A5/3F10 . 1 ,A5) 

J=J+39 

5000 CONTINUE 

ENCODE <5/ 4(?C0 / IBUF<J) )FMT4 
4000 FORMAT <A5) 
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67 ENCODERS, 6000, IBUF<J+1)>FMT3 

70 6000 F0RMAT<A5) 

71 LENGTH=J 

72 CALL STRTFC IPV,NAMER> 

73 CALL APNDFCIBUF<1>, LENGTH) 

74 CALL APNDF<IBUF<J + 1), 1) 

75 CALL CLOSF 

76 ■ RETURN 

77 END 
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2 

3 

4 

5 

6 
7 


SUBROUTINE LOAD<A,T,N) 
DIMENSION T<i)/A<l> 

DO 1 1=1, N 
T<I >=A<I> 

1 CONTINUE 
RETURN 
END 
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1 

2 C 

3 C 

4 C 

5 

6 

7 . 

10 

11 

12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 

46 
4 7 

50 

51 

52 

53 

54 

55 

56 

57 


PROGRAM MAIN 
Ml .0 

MAIN CONTROL FOR IMP 

COMMON/'PARA/ LFLG< 10 > , MFLAG<4> . SUM. LETTER. NO. NFLG. INO. NFLA6C 12) 

* . JJJ.KKK.LLL.MMM.NNN.NAME<2.60).NME<2). IINUM.RRNUM 
COMMON.'GRD/ BBOX < 50 ) . TEST . TEST3 

COMMON/CNT/' I PASS . NUP , NAMER < 3 ) . IPV. ISTRT. I END. IT. LINES 
COMMON^NMES^ NMBB<2) 

COMMON/MOVIT/ NAMDT<2) . XCG. YOG. 2CG. IDTFLG 
COHMON/OPT/ NHIT.OPDPL 

COMMON/'IDB/ NMB<2).IUT<15).TEMP<18). lUJRITE. IJ. ITHMC2) . TIXX< 15) . NUMO 

* . IPRINT.PUIT<3).DCG<60.3). IGNMEC3.60) 

COMMON/'STR/ NTOT . VALUE 1 < 60 ) . VALUE2C60 ) . VALUES C 60 >. VALUE4 < 60 ) . 

* VALUES <60 ) . VALUE6<60 ) . VALUE? <60 ) . VALUES <60 > . VALUES <60 > . 

* VALUEIO <60 ) . VALUE 1 1 <60 > . VALUE12<60 ) .VALUE 13 <60 ) 

REAL VALUE 1 . VALUE2 . VALUES . VALUE4 . VALUES . VALUE6 . VALUE7 . VALUES. 

* VALUES. VALUEIO . VALUE 1 1 . VALUE 12. VALUE 13 
DIMENSION TEST<200) . V<60. 13) . OPDPL < 20 ) . NMRDPL < 2 ) 

EQUIVALENCE <V<1. 1 ). VALUE 1 < 1 ) ) 

INTEGER NAME. NAMER. NMEj IPARA<5) . NMDSPL<2) . IAPN<5) 

LOGICAL IDTFLG 

DATA I INUM/0/. RRNUM/0 / . NME-'2* ' ' MFLAG/4 + 0/'. IDTFLG/ . TRUE . / 

DATA IFARA/IH . 4H UJT . 4HX-CG . 4HY-CG . 4H2-CG/ . I DPV/ 1 1 0 / . I PV/ 11 0 B/ 
DATA IAPN/5H NAME.5H HJT.5H X-CG.5H Y-CG.5H Z-CG/ 

DATA NFLAG/12*0/. LFLG< 10 )/0/. LINES/20/. NMRDPL< 1 )/5HDATA / 

DATA NMRDPL<2)/5H /,P(llT/3*0./ 

DATA2D <BBOX) 

ZSET<0.F) 

MOVE< 1 . F. 1 . F) 

DRA0J< 1 . F. -1 . F) 

• DRAU<-1 .F.-l .F) 

DRAUK-l ,F. 1 .F) 

DRAUKl ,F. 1 .F) 

MOVE< 1 . F. . 7142F) 

DRAW<-1 .F. . 7142F) 

M0VE<-1 .F. ,4285F) 

DRAU/<1 ,F. .42S5F) 

MC7E< 1 ,F. . 142SF) 

DRAUK-l .F. . 1428F) 

M0VE<-1 . F. 1428F) 

DRAUKl .F.-. 142SF) 

MOVE< 1 ,F. - .42S5F) 

DRAU<-1 ,F.-.42S5F) 

M0VE<-1 .F.-.7142F) 

DRAUKl .F.-.7142F) 

M0VE< .7142F.-1 ,F) 
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6 
7 
10 
11 
12 

13 

14 

15 
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31 

32 

33 
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44 

45 
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54 

55 

56 

57 
60 
61 


J.RjM,. «. ,714LT, .7142F) 

MOVc.< . 4285F/ . 7142F) 

. 4285F,-1 .F) 

ACfVE< . H28F,-1 ,F) 

DRAUI C . 1428F. . 7142F) 

MOVEC-. 1428F. . 7142F> 

DRAtii<“. 1428F,“1 ,F> 

M0VE<-.4285F,-1 .F) 

DRAUK- . 4285F, . 7142F) 

M0VE<-.7142F. . 7142F) 

DRAUK-.7142F,-1 .F) 

ENDLIST 

ENODATA 

C 

CLEAR COPDPL) 

30 DISPLAY OPTIONS 
DO 39 1=1,10000 
JS 1XZ=I 
C M.l 

C INITIALIZE PROGRAM VARIABLES 
C 

40 CALL RESET2 
50 CALL RESETl 
100 LABEL <TEST) 

2SET<0 .OF) 

IF<<MFLAG<1) .EQ.O) .AND. <MFLAG<2> .EQ.0>) GO TO 130 
WRITEC 16, 120 ) HMDSPL<1>,NMDSPL<2), IPARAUNO) 

120 FORMAT< "sesX3060aY3056” , 2A5, 2X, A5> 

130 IF'JJJ.NE.O) GO TO 170 

IFCLFLOC 1 ) . EQ , 1 ) GO TO 150 
tlRIT£<16, 140) IINUM 
l'»0- FORMATS ”3B3X.90 70 8Y30 64" ,110) 

GO TO 190 

150 WRITE< 16, 160 ) RRNUM 

160 FORMAT< "383X30 703 YQO 64" , FIO . 5) 

GO TO 190 

170 a)RITE<16, ISO) NME<1),NMEC2) 

180 FORMATC "383X30 70 3V30 64", 2A5) 

190 ENDLIST 

LABELCOPDPL) 

2SET<0 .OF) 

WRITEa6,S10 ) IDPV,NMRDPL<1 ), NMRDPL<2) 

: to FORMAT< " SBsSsXsO 773 YsO 72 " ,14, " S.XSO 773Y30 76 " , 2A5) 
ENDLIST 
C M . 2 

C N:RKhL exit from program 
C 

195 IF< .N0T.DEPRESS<16)> 60 TO 200 
CALL tNHALT 
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3 

4 

5 

6 
7 
10 
11 
12 

13 

14 

15 

16 
17 
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50 

51 

52 

53 
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57 
60 
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62 

63 

64 

65 

66 


CALL STOPS 
CALL EXIT 

200 IF<MrLAG<3) . EQ. 0> GO TO 195 
CALL AN6 
C Ml .3 

C PACK/VOL OPTION 
C 

IF<NFLAG<n .EQ.O) GO TO 500 
MFLAG<2)-1 ' 

NO = l 

flMDSPL<l >=5HENTER 
MMDSPL<2>=4H IPV 
IF<LFLG<3> ,NE, O GO TO 100 
IDPV=I INUM 
IPV=IOCT(I INUtO 
GO TO 40 
C Ml .4 

C FILE NAME OPTION 
C 

500 IF<NFLAG<2) .EQ.O > GO TO 510 
MFLAG<1)=1 
NO=l 

NMOSPL< 1 )=5HENTER 
NMOSPL<2)=5H FILE 
JJJ = 1 

IFCLFLGO) .NE . 1 > GO TO 100 
NAMERCl >=NME< 1 > 

NAMER<2)=NME<2) 

NMRDPL<1 >=NME< 1 ) 

NMRDPL<2)=NME<2> 

GO TO 40 
C Ml . 5 

C LIST VOLUME OPTION 
C 

510 IF<NFLAG<3) .EQ.O) GO TO 520 
NO=l 

HMDSPL<1>=5HTBD 
NMDSPL(2>=5HEHTER 
IF‘'LFLG<3> .NE. 1) GO TO 100 
GO TO 40 
C M 1 . 6 

C BLACK BOX OPTION 
C 

520 IF<NFLAG<4> .EQ,0> GO TO 530 
NFLAGC6>=0 

IF< <MFLAG< 1 ) . EQ . 0 > . AND . <MFLAG<2> .EQ.0>> GO TO 100 

DO 525 1=1.2 4 

NMB<I)=NAME<I.NO> 

525 CONTINUE 

DO 526 J=1 . 13 
UjT<:j> = V<NO. J> 

IF<<J.LT.2> .OR. CJ,GT.4)> GO TO 526 
I=-J-l 

pwr<i)=v<:No. 
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67 

526 

CONTINUE 

70 

527 

iFCLFLG(6.'> .EQ. 1> GO 

71 


CALL (UABCOM 

72 


MAMEC1<N0>=NMBBC1> 

73 


NAME<2>N0)-NMBB<2> 


61 


MAIN 


VER 1 


REV B 


16 JUN 76 


PAGE 4 


1 DO 528 J=1 V 13 

2 V<NO. J)=WT<J) 

3 ir<<J.LT.2> .OR. <J.GT.4>> GO TO 528 

4 I«J-1 

5 DCG<.N0, n^DCCCNO, I) + <V<NO» J>-PWT<I)> 

6 528 CONTINUE 

7 NFLAG<4>=0 

10 NUP=i - 

11 NFLAG<6>=! 

12 GO TO 50 

13 C K1 .7 

14 C GEOMETRY OPTION 

15 C 

16 530 IF<NFLAG<5) .EQ.O) GO TO 540 

17 NFLAG<4)=NFLAG<6>=0 

20 CALL UABCOM 

.21 NT0T = NT0T+1 

22 NAME<1.NT0T>=NMBB<1) 

23 NAME<2, NT0T)=NMBB<2) 

24 DO 535 J=1 , 13 

25 V<NTOT> J>=0IT< J> 

26 535 COHTINUE-- 

27 DO 536 1=1,3 

30 DCGCNTOT, I )=0 . 

31 536 CONTINUE 

32 NFI.AG<5>=0 

33 NFLAG<6)=1 

- ^ 34 NUP=1 

' 35 GO TO 50 

36 C Ml .8 

37 C -DISPLAY OPTION 

40 C 

41 540 IF/NFLAGCG) .EQ .0) GO TO 550 

42 IFCNUP.EQ.l) GO TO 543 

43 IFCIPASS.NE.O) GO TO 544 

44 NAMER<1>=NMRDPL<1) 

45 NAMERC2)=NMRDPL<2) 

46 CALL LOCFLECIPV.NAMER, IFILNO) 

47 OPEN<20 , IPV, IFILNO) 

50 READ <20, 541) NTOT 

51 541 F0RMAT<I3) 

52 READ<20,542) <NAME(1, I),NAME<2, I ) , VALUEl < I ) , VALUE2< I ) , VALUE3< I ) , 

53 ♦ VALUE4< I ) , VALUE3< I ) , VALUE6< I ) , VALUE7< I > , VALUESC I ) , VALUE9< I ) , 

54 ♦ VALUE10<I),VALUE11<I),VALUE12<I),VALUE13<I), IGNME<1,I), 

•55 * IGNME<2, I), <DCG<I,J),J=1,3),I=1,NT0T) 

56 542 FORMAT<2A5,6F10 . 1,>',7F10 . 1,,^, 2A5,3F10 . 1) 

57 CLOSE<20> 

60 543 CALL WABSET 

61 544 ITEMP=NTOT/LINES 

62 IRMDR=NTOT//LINES 

63 IF<IRMDR.NE.O) ITEMP=ITEMP+1 

64 IF< .NOT .DEPRESS<32>) GO TO 546 

65 IT=IT+1 
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IP«SS=0 

5^6 IFUT.GE. ITEMP> IT=0 
ISTRT=1+(IT*LINES) 

IEND=LINES+<IT*LINES> 
irCtEND.CT.NTOT) IEND=NTOT 
CALL UJABLST 

IF<<MFLAG<D ,EQ. n .OR. CHFLAG<2> .EQ. 1>) CO TO 547 
DO 545 I = D 12 
IF<I,EQ.6) GO TO 545 
IF<NFLAG<D .NE. D GO TO 545 
NFLAG<65=s0 
GO TO 50 
545 CONTINUE 
547 NMDSPL<D=NAME<DNO) 

N?'1DSPL<2>=NAME<2.N0) 

IF<LFLG<3) .EQ.0> GO TO 100 
GO TO 50 
Ml .9 

C APPEND DISPLAY OPTION 
C 

550 IF<NFLAG<7> , EQ. 0> GO TO 560 
IF<IX.EQ.6) GO TO 559 
MFLAG<2)=1 

IF<IX.NE. D GO TO 552 
MFLAG< D=JJJ=1 
552 NMDSPL< 1 >=5HENTER 
NHDSPL<2)=IAPN< IX) 

IF<LFLG<3) .NE. 1 ) GO TO 100 
IF<IX.NE. 1) GO TO 555 
NTOT=NTOT+l 
NAME<1>NT0T)=NME<1) 

' HAME<2, NT0T)=NME<2) 

GO TO 557 
555 J=IX-1 

V<NTOT, J)=RRNUM 
IF<J .EQ. 1 ) GO TO 557 
I = J-1 

DCG<NTOT. I)=0 . 

557 IX=IX+1 
GO TO 50 

559 DO 55S J=5, 13 
V<NT07. J)=0 . 

555 CONTINUE 

DO 556 I = D2 
IGNME< D NT0T)-5H 

556 CONTINUE 
IPASS=0 
NUP=1 
IX=1 

NFLA6<7>=0 


m 
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• 

1 


NFLAG<6)=1 


2 


CO TO 50 


3 

C 

Ml . 10 


4 

C 

DELETE COMPONENT OPTION 


5 

c 



6 


560 IF<NFLAG<8) .EQ.0> GO TO 

570 

7 


IF<<MFLAG<1 > .EQ.O) . AND. ' 

CMFLAG<2> .EQ 

10 


NT0T=NT0T-1 


11 


DO 566 I=N0.NT0T 


12 


DO 564 J=l,2 


13 


NAME<Jy I)=NAME<J. I+l> 


14 


IGNME<J. n = IGNME<J, I + l) 


15 


564 COrjTINUE 


16 


DO 565 K=1 , 13 


17 


V<I.IO=V<I + l,K) 


20 


IF<<K.LT.2> .OR. <K.GT.4)5 

' GO TO 565 

21 


N=K-1 


22 


DCG<I/H)=DCG<I+1>N) 


23 


565 CONTINUE 


24 


566 CONTINUE 


25 


IPASS=0 


26 


NUP=1 


27 


NFLAG<8>=0 


30 


NFLAG<6)=1 


31 


GO TO 50 


32 

z 

M. . 1 1 


33 

c 

stokl option 


34 

c 



' 35 


570 IF<NFLAG<9> .EQ.O) GO TO 

580 

36 


NAMER< 1 )=NMRDPL< 1 ) 


37 


NAMER<2)=NMRDPL<2> 


40 


-CALL ST0RE< IPV.NAMER) 


41 


NFLAG<9)=0 


42 


HFLAG<6)=1 


43 


GO TO 50 


44 

c 

Ml .12 


45 

c 

INPUT NEW DATA FILE OPTION 


46 

c 



47 


580 IF<NFLAG<10> .EQ.O) GO TO 

590 

50 


CALL INPUT<IPV) 


51 


GO TO 30 


52 

c 

Ml . 13 


53 

c 

BLACK BOX TRACKER OPTION 


54 

c 



55 


590 IF<NFLAG<11) .EQ.O) GO TO 

595 

56 


IF<<MFLAG< 1) .EQ.O) .AND. < MFLAG < 2 > . EQ . i 

57 


NAMDT<l)=NAME<l/HO) 


60 


NAMDT<2)=NAME<2.N0) 


61 


XC6=V<N0.2) 


62 


YCG=V<N0.3) 


63 


2CG=V<N0.4) 


64 


DO 593 J=2.4 


65 


I = J-1 


66 


PUIT<I)=V<NO. J) 
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GO TO 100 
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67 593 CONTINUE 

70 CALL CGMOVE 

71 V'<N0.2)=XCC 

72 VCN0*3)=VCG 

73 VCN0>4>*ZCG 
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5 

6 
7 
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12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 


37 

40 

41 

42 
45 

44 

45 

46 

47 
50 


DO 594 J=2,4 
I*J-1 

DCG<NO, I )=DCGCNO, I )+< V<NOi J>-PUT< I ) ) 

594 CONTINUE 
NFLAGai)=0 
NUP=1 

NFLAG<6>=1 
GO TO 50 
C Ml. 14 
C PRINT OPTION 
C 

595 IF<NFLAG<12) .EQ.O) GO TO 100 
CAlL iPPLCOPY 

NFLAGC 12)=0 
NFLAG<6)=1 
GO TO 50 

SUBROUTINE RESETl 
C M2 .0 

C FIPsT INITIALIZATION OF MAIN PROGRAM VARIABLES 
C 

MFLAG< 1 )=MFLAG<2>=LFLGC10 )=NO=NFLG=JJJ=0 

LLL=KKK = MMM=NNN=IUlRITE=IPRINT = 0 

INO=l 

NUM0=9 

DO 600 1=1,3 
600 PUIT<I> = 0. 

RETURN 

SUBROUTINE RESETS 
C H3.0 

C SECOND INITIALIZATION OF MAIN PROGRAM VARIABLES 
C 

. MFLAG<3) = I INUM=IPASS=IT=0 
IX=1 

RRNUM = 0 , 

NMDSPLO ) = NMDSPL<2)=5H 
DO t I = 1 , 1 2 
NFLAG<I>=0 
1 CONTINUE 
RETURN 
END 
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1 SUBROUTIME MXM3X3<A. B/ C) 

2 C 

3 C C=AB ALL 3X3 MATRICES 

4 C 

5 DIMENSION A<9> > B<9) , C<9> 

6 C< 1 >=A< 1 )*B<: 1 )+A<4>^BC2)+A<7)*«'B<:3) 

7 C<2)=AC2>-*<B< 1 > + A<5)*B<2)+A<8)+B<3) 

10 C<3)=A<3>*B< 1 ) + A<6)=+=8<2)+A<9)*B<3> 

n C<4)=A< 1 )‘*«B<4)+A<4)*B<5)+A<7>*B<6> 

12 C<5>=A<2)+B<4)+A<5)*B<5)+A<8:>*B<e> 

13 C<6)=A':3> + B<4)+A<6>*B<5)+A<9>>(‘B<6) 

14 C<7)=A<1 >-+'B<7>+A<4)*^B<8)+A<7)+B<9) 

15 C<8)=A<2)*B<7) + A<5)‘4'8<8>+A<8> + B<9> 

16 C<9)=A<3>*B<7)+A<6)*B<8)+A<9>*B<9> 

17 RETURN 

20 END 
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1 SUBROUTINE MXV<M,V.VV> t :• 

2 C 

3 C THIS ROUTINE MULTIPLIES A 3X3 MATRIX TIMES A 3Xl MATRIX AND RETURNS 

4 C THE RESULTING 3X1 MATRIX. 

5 C 

6 C INPUTS! M = 3X3 INPUT MATRIX 

7 C V « 3X1 INPUT MATRIX <VECTOR> , ! 


10 ' 

C 



11 

C 

OUTPUT! VV = 3X1 

OUTPUT. MAY BE STORED IN ’SAME LOCATION 

12 

C 

AS 

INPUT VECTOR. i. j 

13 

C 


• • 1 

14 


DIMENSION V<3>.TMP<3) 

.VVC3) * * 1 

15 


REAL M<9) 

1 

1 

16 


TMP<I> = M<n*v<i) + 

MC4)*VC2) + MC7)'*'VC3) ' ’ ’ : 

17 


TMPC2) = M<:2)>«‘V<1> + 

MC5>*VC2) + MC8)«VC3) • ' 

20 


TMP<3) = M<3)*V<n + 

M<6)*VC2) + MC9)*V<3) 

21 


DO 5 1=1.3 

1 

22 


VV<I> = TMPC I? 


23 


5 CONTINUE 


24 


RETURN 

• 

25 


END 
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SUBROUTINE PAT<A,U.R,B) 

C 

C TRANSFORMS INERTIAS BY PARALLEL AXIS THEOREM 
C 

DIMENSION AC6)/R(3>.B<6? 

B< 1 >=A<1 )+iJJ*<R<2)**2+R<3)'*'*2> 
B<2)=A<:2)+liJ*(:R< 1 )*>*<2+R<3)**2) 
B<3>=A<3)+U1h«<R< 1 >+*2 + R<:2>*‘*<2> 

B<4 5=A(4)+UJ>t<R<2)>*<R< 1 > 

B<5)=A<5> + U*R< l)t=R<3> 

B<6> = A<6)+U*R<:3)*RC2) 

RETURN 

END 
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22 
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26 
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30 

31 
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33 

34 

35 

36 

37 
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41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 
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F INCLUDE=FORMT 

SUBROUTINE STORE < IPV, NAMER> 

C SI .0 

C STORE DATA TO DISK 
C 

COMMON/PARA/' LFLGC 10 ) . MFLAGC4) > SUM. LETTER y NO. NFLG, INO. NFLAG< 12) 

♦ / JJJ. KKK.LLL. MMM.NNN. NAME <2. 60 ) . NME<2> . I INUM. RRNUM 
COMMOH/'STR/ NTOT . VALUE 1< 60 > . VALUE2<60 > / VALUES <60 ) . VALUE4<60) , 

♦ VALUES < 60 ) . VALUE6C60 ) . VALUE7< 60 >. VALUES < 60 > . VALUES <60 > . 

♦ VALUE 10 <60 ) . VALUEl 1 <60 ) . VALUE 12<60 ) . VALUE13<60 ) 

COMMON/'tllB/' NMB<2) . UT< 15> . TEMP< 18? . lURITE. IJ. I THM < 2 > . T I XX < 1 5 > . NUMO 

♦ . IPRINT. P1UT<3) . DCG<60 . 3) . IGNME<3. 60 ) 

DIMENSION IBUFOOOO? 

REAL VALUEl. VALUE2. VALUES . VALUE4 . VALUES . VALUES , VALUE? . VALUES . 

♦ VALUES. VALUEl 0 . VALUEl 1 . VALUE 1 2 . VALUE 1 3 
INTEGER NAME. NAMER<3> 

DATA FMT2/-1515151515B/ / FMT3/7777777777B/'. FMT4/'7777700000B/' 

CALL BLDNAM<NAMER) 

ENCODE<5. 1000 . IBUF< 1) ) NT0T.FMT2 
1000 F0RMAT<I3. A2) 

J=2 

DO 5000 I=l.NTOT 

ENCODE<75. 2000 . IBUF<J)> NAME<1. I),NAME<2. I ) . VALUEl < I ? . VALUE2< I ) . 

♦ VALUE3< I ? . VALUE4< I ) . VALUE5< I ) . VALUE6< I ) . FMT2 
2000 F0RMAT<2A5. 6F10 . 1 . AS) 

ENCODE <75. 3000 . I BUF < J+ 1 5 ) ) V ALUE7 < I ) . VALUES < I ) . VALUES < I ) . 

♦ VALUEIO < I ) . VALUEl 1 < I ) . VALUE12< I> . VALUEl 3< I ) . FMT2 
3000 FORMAT<7F10 . 1 . A5) 

ENCODE<45. 7000 . IBUF< J+30 ) ) IGNME<1 . I ) . IGNME<2. I). <DCG< I .K). 

♦ K=1 . 3>. FMT2 

7000 F0RMAT<2A5. 3F10 . 1 . AS) 

- J=J+39 
5000 CONTINUE 

ENCODE <5. 4000 . IBUF< J) )FMT4 
4000 F0RMAT<A5) 

ENCODE <5. 6000 . I BUF < J+ 1 ) ) FMT3 
6000 F0RMAT<A5) 

LENGTH=J 

CALL STRTF< IPV. NAMER) 

CALL APNDF< I BUF <1). LENGTH) 

CALL AFNDF< IBUF< J + 1 ) , 1 > 

CALL CLOSE 

RETURN 

END 
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SUBROUTINE TOKAiB.N) 
DIMENSION A<n,BU> 

DO 1 1=1. N 
ft<I)=ft<I)+B<n 
B< I)=0 . 

.1 CONTINUE 
RETURN 
END 
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5 
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25 

26 
27 

30 

31 

32 

33 
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35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

55 

56 

57 
60 
61 
62 

63 

64 

65 

66 


r INCLUDE=DISK.TTY,PPL 

5UBr:QUTINE UIABCOM 

c mi.o 

C MASS PROPERTIES CALCULATIONS 
C 

COMMON/'PARA/ LFLG< 10 MFLAG<4) , SUM> LETTER. NO. NFLG, IHO. NFLAG< 12) . 

* JJJ.KKK.LLL.MMM.NNN.NAME<2.60).NME':2). IINUM.RRNUM 
COMMON/UB/' NMB<2) . UJT< 15) . TEMP< 18) . lURITE. IJ. ITNM(2) . TIXXC 15) . NUMO 

* . IPRINT.PUJT<3;.DCG<60.3). IGNMEC3.60) 

COMMON/'GRD/ BBOXCSO ) . TEST<200) . TEST3<50 ) . OPDPL<100 ) . BOX <20 ) 
COMMON/LTR/' IFLAG<40) 

COMMON/'NMES.' NMBB<2) 

COMMON/STR/ NTOT . VALUE 1 < 60 ) . VALUE2<60 ) . VALUES <60 ) . VALUE4<60 ) . 

* VALUES <60 > . VALUE6<60) .VALUE? <60 ) .VALUES <60 ) . VALUES < 60 ). 

* VALUEIO <60 ).VALUE1 K60). VALUE12<60 ) . VALUE 1 3 < 60 ) 

D IMENS ION L 1 < 3 ) . X I N < 4 ) . Y I N < 4 ) . Z I N < 4 ) . L2 < 3 ) . D 1 < 3 ) . 2U < 3 ) 

DIMENSION YU<3) .XU<3) , TEMP2<9) . DC < 9 ) , PM I P < 9 ) . PM I < 9 ) 

DIMENSION DELX<3) . XO <3) . XBAR<3) . XBRF < 3 ) . G I N < 1 0 0 ) 

DIMENSION SIXX< 15) . D2R<9> . XA < 30 ) . B I N < 1 50 ) 

DIMENSION YA<30).ZA<30).XB<30). YB<30).ZB<30) 

DIMENSION FACT<10).NMDS<6). IFILE<3). IFILENO) 

EQUIVALENCE < I YY . UT < 6 ) ) . < I ZZ . UIT < 7 > ) . < I XY . UJT < 8 > ) . < I XZ . WT < 9 ) ) 
EQUIVALENCE < I YZ . UIT < 1 0 ) ) . < S . UIT < 1 1 ) ) . < IXX.UT<5)) 

EQUIVALENCE < AREA . WT < 1 3 ) ) . < VOL . lit)T< 1 2) ) . < DELY , DELX < 2 ) ) 

* . <DEL2.DELX<3)) 

EQUIVALENCE <FACL,FACT< 1 ) > . < FAC I . FACT < 2 ) ) . < FACW . FACT < 3 ) ) 
EQUIVALENCE < FACG . FACT < 4 ) ) . < PS I . FACT < S > > , < THETA . FACT < 9 ) ) 
EQUIVALENCE < PH I . FACT < 1 0 ) ) 

REAL LI . L2. IXX. lYY. IZZ. IXY. IXZ. lYZ 
LOGICAL IFLAG.SUB.BLKBX.RFLAG.AFLAG.BFLAG 
INTEGER STAT. ROOT. STNAME.STATT. OVL 
C 

DATA DELX/-3*0 ./•PSI.PHI.THETA/'3+0 ./RHO.H/S*! .x 

DATA FAC I /I ./.FACUJ/'l ./.FACG/l . / . FACL/ 1 2 . z' . TEMP/ 1 8 + 0 ./ 

DATA IFLAG/40 + . FALSE. /.TlXK/15+0 IPRINT/0/. IBASE/0/ 

DATA NMB< 1 )/4HC0MP/. STNAME/4HSUBT/ . I WR I TE/0 / . IREFL/O/ 

DATA tUT/15+0 SUB/. FALSE. /. IFLP/0/ 

DATA NME/’ '/.NMB<2)/' ' /. OVL/1 lOB/. I IPV/1 108/ 

DATA NMDS/5HENTER. 5H VALU.5HE .5HENTER.5H NAME.5H / 

DATA IFILE/5HDATA .5H . 5H / 

C 

C FORMAT STATEMENTS 
C 

1000 FORMAT</. 1X.2A5./6H UT = .G13.7./6H XCG =.G13.7.5X. 

* 5HYCG =. G13 . 7. 5X. 5HZCG =.G13.7./6H IXX =.G13.7.5X. 

1- 5HIYY =. G13 . 7. 5X. 5HIZZ =.G13.7./6H IXY =.G13.7.5X. 

+ 5HIXZ =.G13 .7.5X.5HIYZ =. G13.7) 

1011 F0RMAT<6H ASF = . G 1 3 . 7 . 5X . 5HV0L =. G 1 3 . 7. 5X. 5HAFR =.G13,7) 

150 0 FORMAT<2<3F10 ,0. ID) 

1700 FORMAT <A5) 

1710 F0RMAT<I5) 

C 

C - 

NUM0 = 9 
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67 IF<IJ.CT.1> GO TO 6000 

70 NTIME=0 

71 6000 IFC ,N0T.SWITCH<1>> GO TO 6660 

72 IF<IIURITE.EQ. 1> CLOSE<20) 
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31 

32 

33 
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C U1 . 1 

C INITIALIZE SUBROUTINE VARIABLES 
C 

6660 CALL RESET 

IF<NFLAG<6> .NE. 1) GO TO 4 
CALL RESET69 
BLKBX=.TRUE, 

GO TO 77 

4 IF<NFLAG<4) .NE. n GO TO 5 
SUB= , TRUE. 

BLKBX=.TRUE. 

IFLP=1 
GO TO 78 
C Ull .2 

C GEOMETRY INPUTS SETUP 
C 

5 BLKBX=. FALSE. 

DISPLAY GEO 
CALL RESETS? 

6 CALL RESET 
LABELCGIN) 

ZSET<0 .OF) 

WRITE < 16. 109) NME<1).NMB<2>. IDPV. IFILECl), IFILE<2>.RH0 

108 FORMATC "SBSXa0 50SiYs0 30 " .2A5. 

* "8XSi050aYa034 " . 14. 

* "OXa050sY3040" . 2A5. 

* '‘9X80503Y9044'' . FIO , 3) 

WRITEaS.lOS) H. FACL. FACI . FACW 

109 FORMAT<"3BaX3050aY8050‘'.F10 .3. 

* "8X30508Y3054" . no . 3. 

* "0X3O5O8Y9O6O " . FiO . 3. 

* "3X30503Y3064" . FIO . 3) 

WR1TE<16. 1 10) FACG. WTF. IREFL. I WRITE. SUB 
no FORMAT <" 808X30 50 3 Y30 70 " . F 1 0 . 3 . 

* "3X30503Y3074" .FIO .3. 

* "8X31403Y3030 ''.-II. 

* "3X31403Y3034 ".II, 

* "3X31403Y3040 " . L9) 

WRITEC 16. 1 1 1 ) DELX< D.DELY.DELZ.PS I 

1 1 1 F0RMATC"3B3X31403Y3044" . FIO ,3. 

* "3X31403Y3050 " .FIO .3. 

* "3X31403Y3054" . FIO . 3. 

* "3X31403Y3060 " .FIO .3) 

WRITE< 16. 1 12) THETA. PHI . IBASE 

112 r0RMATC"3B3X31403Y3064".F10 .3. 

* "8X.31403Y.30 70 ".FIO .3. 

* "8X31408Y8O74 ">ID 

ENDLIST 

15 LABELCTESTS) 

ZSET<0.0F) 

IFC .NOT. IFLAGOO) ) GO TO 99 
1 = 1 

IFCIFLAGOD) 1=4 

WRITE<16.92) NMDS<I).NMDS<I+l).NMDS<I+2) 
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67 

70 

71 

72 

73 

74 

75 

76 

77 
100 
101 
102 

103 

104 

105 

106 
107 

no 

111 

112 

113 


92 rORMAT<"OBsXS0609Y9056", 3A5> 

89 lF<IFLAG<3n> GO TO 96 
IF<LrLG<l) .EQ. 1) GO TO 83 
UIRITE<16,90> IINUM 

90 FORMAT < "889X60 70 eVeO 64". 110) 

GO TO 93 

88 WRITEC16. 91) RRNUM 

91 FORMAT<"8BaX90709Y9064".F10 .5) 
GO TO 93 

96 WRITE <16. 97) NME < 1 ) . NME < 2) 

97 FORMAT< "8B9X90708Y9064". 2A5) 

93 ENOLIST 

IF< .NOT. IFLAGOO)) GO TO 15 
CALL AN7 

IF< . NOT. IFLAG<3) ) GO TO 320 
IF<LFLG<3) ,NE. 1 ) 60 TO 15 
ENCODE <5.. 170 0 .NMB<D) NMECl) 
ENCODE <5. 1700 .NMB<2)) NME<2) 

GO TO 6 

320 IF< . NOT. IFLAG<4) ) GO TO 321 
IF<LFLG<3) .NE. 1) GO TO 15 
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IDPV=IINUM 

IIPV*IOCT<IINUM> 

GO TO 6 

i4.1 IFC .NOT. IFLAGC5)) GO TO 322 
IFCLFLGO) ,NE. 1) GO TO 15 
ENCOOE<5/ 1700, IFILEd)) NME<1> 
ENC0DEC5, 1700> IFILE<2>) NME<2) 
IGNMEa>NTOT + l> = IFILE<n 
IGNME<2.NT0T+1)=IFILE<2> 

00 3200 1=1.3 
3200 IFILEN«:i) = IFILE<I) 

CALL LOCFLECI IPV. IFILEN. IFILNO) 
GO TO 6 

322 IF< . NOT. IFLAG(:6>) GO TO 325 
IF<LFLG<3> . NE . 1) GO TO 15 
RHO=RRNUM 
GO TO 6 

325 1F<.N0T. IFLAG<7>) GO TO 326 
IF<LFLG(3> . NE . 1 ) GO TO 15 
H=RRNUM 

GO TO 6 

326 DO 327 1=1.4 
J=I + 7 

IF< .NOT. IFLAG<J>> GO TO 327 
IFCLFLGCS) .NE. 1) GO TO 15 
FACT<n = RRNUM 
GO TO 6 

327 CONTINUE 

IF< .NOT. IFLAGC12?) GO TO 328 
IFCLFLGCS) .NE. 1) GO TO 15 
IREFL=IINUM 
■ GO TO 6 

328 IFC .NOT, IFLAG<13>> GO TO 329 
IF<LFLG<3> .NE, U GO TO 15 
IWRITE=I INUM 

GO TO 6 

329 IFC . NOT. IFLAG-:i4>) GO TO 330 
IF<LFLG<3> . NE . 1> GO TO 15 
SUB=. FALSE. 

IF<NME<n . EQ . IHT) SUB=.TRUE. 

• GO TO 6 

330 DO 331 1=1.3 
J=I+14 

IF< .NOT. IFLAG<J>) GO TO 331 
IF<LFLG<3> .NE. 1) GO TO 15 
DELX<n = RRNUM 
GO TO 6 

331 CONTINUE 

DO 332 1=8. 10 
J=I+10 

IF< . NOT. IFLAGC J) ) GO TO 332 
IF<LFL6<3) .NE. 1) GO TO 15 
FACT<I>=RRNUM 
GO TO 6 
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67 332 CONTINUE 

70 IFC.N0T.irLAG<22>> GO TO 359 

71 ir<LFLG<3> .NE. 1) GO TO 15 

72 6JTF=RRNUM 

73 GO TO 6 

74 359 IF< .NOT. IFLAG<23>> CO TO 358 

75. IF<LFLG<3) .NE. 1) GO TO 15 

76 IBASE = IIHUfi 

77 GO TO 6 

100 358 IF<IFLAG<2D) GO TO 77 

101 GO TO 15 

102 c un .3 

103 C BLACK BOX INPUTS SETUP 

104 C 

105 360 IFLP=0 

106 CALL RESET69 

107 IFILE<l)=IGNME<l.NO) 

110 IFILE<2)=IGNME(2/N0) 

111 DISPLAY BKBX 

112 797 CALL RESET 

113 LABEL<BIN) 

114- • • ZSETCO .OF) - 

115 UJRITE< 16.208) NMB< 1 ) . NMB<2) . <UIT< I ) > 1 = 1 . 7) 

116 208 FORMATC“6iB3X3030aY3030".2A5. 

■ 117 ♦ '•QXa030QY9034“.F10 .0, 

•' 120 * ••aXQ030aY3040".F10 .0. 
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7 
10 
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21 

• 22 

23 

24 

25 

26 
27 
30 

. 31 
.* 32 

33 

34 
*35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 
, 55 

56 

57 
60 
61 
62 
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64 

65 

66 


* "SXSl 1O8Y0O4O •* , FIO . 0 * 

* «8XQ170qYQ040",F10.0/ 

* "OX8030QYQ044" ,F10 . 0. 

* ''8X81108Y8044'SF10 ,0> 

* "8X81708Y8044",F10 . 0) 

tURITE<16.212> 1 = 8. 13> 

212 FORMAT <"8B9X80 308Y8050'*.F10 .0. 

* ''8XS110SY8050">F10 .0. 

* "8X8170SY8050".F10 .0. 

* "8X80308Y8054".F10 .0. 

* "aXall08Y8054".FlO .0. 

* "8X81708Y8054''.F10 ,0> 

WRITE<16.209) FACL . DELX < 1 ) . FACI . DELY 

209 FORMAT < "868X8050 8Y80 64 ">F10 .0. 

* "8X81408Y8064" . FIO . 0 . 

* ‘•8X8050QY8070 ".F10 ,0. 

* "OX81408Y8070 " . FIO . 0 J 

yjRITE< 16. 210 ) FACUJ.DELZ. FACG.PSI 

210 FORMAT<"8B8X80508Y8074".F10 ,0. 

* "8X8i40aY8074''.F10 .0. 

* "OX80508Y8100 " . FIO . 0 . 

* "axsMOaYaioo " . Fio . o > 

lURITE<16.2in OVL. THETA. I WR I TE . PH I . I PR I NT . I BASE. IFILE< 1 ) . IFILE<2> 

211 FORMAT<"8B8X8050SYal04 ".03. 

* "8X3l40aY8l04".FlO .0. 

* ''SX80508YO1 10 ".II. 

* "8X8l40SY8liO".F10.0. 

* "«i'X80 50 8Y3l 14 ".II. 

"SX81408Y81 14 “>I1. 

* "8X80508Yal20".2A5) 

ENDLIST 

798 ,LABEL<TEST3) 

2SETC0.0F> 

1F< .NOT. IFLAG<30>) GO TO 213 
1 = 1 

IF< IFLAGOl > ) 1 = 4 

URITEC 16. 214) NMDS < I > . HMDS < I + 1 ) , NMDS < I +2 ) 

214 FORMAT<"8B8X80608Y8056". 3A5> 

213 IF< IFLAG<31 ) ) GO TO 220 

IF<LFLG<1).EQ.1> GO TO 215 
litfRITE< 16. 216) 1 1 HUM 

216 FORMAT<"8B8X80708Y8064" . 110) 

GO TO 218 

215 IJJRITEC 16. 21?) RRHUM 

217 FORMAT<"8E8X80?08Y8064".F10 .5) 

60 TO 218 

220 URIT£< 16. 221 ) NME < 1 ) . MME < 2 ) 

221 FORMAT <"8B8.X80 70 8Y8C 64 ".2A5) 

218 ENDLIST 

219 IFC .NOT. IFLAGOO)) GO TO 219 
CALL AN7 

IF< .NOT. IFLAG<29)) GO TO 799 
1F<LFLG(3) . NE . 1 ) GO TO 798 

■ ENC0DE<5. 1700.NMBC1)) NMECl) ^ 

m 


78 
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67 

70 

71 

72 

73 

74 

75 

76 

77 
100 


ENCG0E<5, 1700.NMB<2>) NME<2) 
CO TO 797 

799 DO 800 1 = 1 / 13 
J=I+2 

IF< .NOT. IFLAG<J>) GO TO 800 
CLFLG<3) .NE. 1) CO TO 793 
UITCI )=RRNUM 
CO TO 797 

800 CONTINUE 

DO 805 1=1.4 
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23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 
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53 

54 

55 

56 

57 
60 
61 
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63 

64 

65 

66 


IF< .NOT. IFLAG(J>) GO TO 805 
IF<LFLG<3) .NE. 1 > GO TO 798 
FftCT<n = RRNUM 
GO TO 797 
805 CONTINUE 

IF< .NOT, IFLAGC20 >) GO TO 810 
IF<LFLG<3> . NE. 1) GO TO 798 
SUB=. FALSE, 

IFCNME<1) ,EQ. IHT) SUB=.TRUE. 

GO TO 797 

810 IF< .NOT. IFLAG<21>> GO TO 815 
IFCLFLGO) . NE. 1 > GO TO 798 
ItURITE=I INUM 

IF<I1URITE.EQ, 1) OPEN<20>OVL) 

GO TO 797 
815 00 820 1=1.3 
J=I+21 

1F< .NOT. IFLAGUJ) GO TO 820 
IF<LFLG<3) . NE . 1 > GO TO 793 
DELXC I >=RRNUM 
GO TO 797 
820 CONTINUE 

DO 825 1=8.10 
J=I+17 

IFC .NOT. IFLAG<J5> GO TO 825 
IF<LFLG<3> .NE. 1> GO TO 793 
FACTCn=RRNUM 
GO TO 797 

825 CONTINUE 

IF< .NOT. IFLAG<32)> GO TO 826 
IF<LFLG<3) .-NE-r-n GO TO 798 
- IPRINT=IINUM 
NUM0=9 

IFCIPRINT.EQ. 1) NUM0=13 
60 TO 797 

826 IF< .NOT. IFLAG<33>) GO TO 823 
IFCLFLG<3) .NE. n GO TO 798 
IBASE=IINUM 

GO TO 797 

828 IF< . NOT. IFLAG<34) ) GO TO 827 
IF<LFLG<3> . NE . n GO TO 798 
ENC00E<5. 1700. IFILEC 1 > > NMEU) 
ENCODE <5. 1700 . IFILECS) > NMEC2) 
ENCODE C5. 1700. IGNME<l.NO>> NME<1) 
ENCODE C 5. 1700. IGNME<2.N0>> NME<2> 
GO TO 797 

827 IF< IFLAG<28> > GO TO 77 
GO TO 798 

C Wl .4 ' 

C COMPONENT COUNTER 
C 

77 NTIME=NTIME+1 
.78 IF<NUM0.EQ.9> CALL ^PTCH9 
C W1 . 5 
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67 C /EJ.0 ARRAYS 

70 C 

71 DO 6001 1=1, 15 

72 TEMPCI)=0. 

73 SIXX<I>=0. 

74 6001 CONTINUE 
75, C tUl,6 

76 C SET UP DC MATRIX FOR DATA TO REFERENCE COORD SYSTEM 

77 C 

100 CALL XFORMCPSI, THETA, PHI, D2R) 

101 C 
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IF<BLKBX) GO TO 7000 
C 

C W1 . 7 

C READ IN ALL SURFACE DATA 
CALL 4NHALT 
NUMI=20 

OPEN<20, IIPV. IFILNO) 

C READ BLANK CARD IN GEOMETRY FILE 
READCNUMI. 1550 ) IDUMMY 
1550 FORMATCIl) 

C 

READ tNUMI. 1500.END=400) X>Y,2,STAT. XX. YY, 2Z. STATT 
RFLAG = .FALSE. 

GO TO 80 

30 IF <RFLAG) GO TO 50 


82 


UABCOM 


VER I 


REV A 


16 JUN 76 


PAGE 7 


1 


RFL4G = .TRUE. 



2 


K * XX 



3 


Y = YY 



4 


2 = 22 



5 


STAT => STATT 



6 


GO TO 60 


• 

7 

50 

. RFLAG = .FALSE. 



10 


READ CNUMI, 1500,END=400> X,Y>2.STAT> 

XX. YY.i 

11 

60 

IF (STAT .EQ. 0 .OR. STAT . EQ . 3) 

GO 

TO 180 

12 


IF (STAT .EQ, 2> GO TO 200 



13 


IF (.NOT. AFLAC) GO TO 200 



14 


MC = M 



15 

80 

M = 1 



16 


IF (STAT .EQ. 2) GO TO 150 



17 


IF (.NOT. BFLAG) GO TO 84 



20 

75 

DO 81 J =1>MC 



21 


XA(J) = XB(J> 



22 


YA(J) = YB(J) 



23 

81 

2A(J) = 2B(J) 



24 

83 

XB(0 = X 



25 


YB(1) = Y 



26 


2B(1) = 2 



27 


GO TO 30 



30 

84 

IF (AFLAC) GO TO 85 



31 


BFLAG = .TRUE. 



32 


GO TO 75 



33 

35 

AFLAC = .FALSE. 



34 


GO TO 83 



35 

150 

AFLAC = .TRUE. 



36 


BFLAG = .FALSE. 



37 


N = N+1 



40 

160 

-XA(M) = X 



41 


YA(M) = Y 



42 


2A(M) = 2 



43 


GO TO 30 



44 

180 

M = M + 1 



45 


IF (AFLAG) GO TO 160 



46 


XB(M) = X 



47 


YB(M) = Y 



50 


2B(M) = 2 



51 


IF (STAT .NE. 3) GO TO 30 



52 

20 0 

MMIN = MINO (M.MC) - 1 



53 


MC = M 



54 

C lUl 

8 



55 

C BE 

GIN COMPUTATION OF SURFACE ELEMENT 

CHARACTERI 

56 

C 




57 


DO 2000 J= l,MMIN 



60 


IPASS=0 



61 


XIN(l) =+DELX( 1 )+XA(J ) 



62 


XIN(2) =+DELX( 1 )+XA(J +1) 



63 


XIN(3) =+DELX(l )+XB< J +1) 



64 


XIN(4) =+OELX(l )+XB( J ) 



65 


YIN(l) =+DELY+YA(J ) 



66 


YIN(2) =+DELY+YA(J +1) 
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VINO) = + DELV + YB<J +1) 

YIN<4) =+DELY+YB<J > 

ZIN<1) =+DELZ+ZA<J ) j 

ZINC2) =+0ELZ+ZA<J +1) s 
ZINO) =+DELZ + ZBCJ +1) 

ZIN<‘4) =+DELZ + ZB<J ) 

TRANSFORM DATA TO REFERENCE COORD SYSTEM 

L1<1)=XIN<‘4)-XIN<1) 

L1<2')=YIN<4)-YIN(1) 

LK3)=ZIN<4>-ZIN<1> 

' • i 

L2<l>=XIN<2)“XIN<i> - 

L2<2)=YIN<2)-YIN<1 ) f 
L20> = ZIN<2)-?I^J< 1 ) 

Dl< 1 )=XINc:4)-XIN<2) 

DK2)=YIN<4)-YIN<2) 

D1 C3) = ZIN<4)-2IN<2) 

XO < 1 )=XINC 1 ) 

X0<2>=YIN<1> 

XO 0)=ZIN< 1 ) 

GO TO 5010 


33 

5020 

CONTINUE 

34 


LKl )=XIN<2)-XIN<3) 

- 35 


LK2)=YIN':2)-YIN<3) 

36 


L1<3)=ZIN<2>-ZIH<3) 

37 

C 


40 


L2< 1 >=XIN<4)-XirK3> 

41 


L2<2)=YINC4)-YIN<3) 

42 


L2C3>=2IN<4)-ZIN<3) 

43 

C 


44 


DK 1 >=XINC2)-XIN<4) 

45 


D1 C2)=YIN<2)-YIN<4) 

46 


D1 ^3)=2IN<2)-ZIN(4) 

47 

c 


50 


X0< 1 )=XIN<3) 

51 


X0<2)=YIW<3) 

52 


X0<3)<;::ZIN<3) 

53 


GO TO 5010 

54 

c 


55 

5010 

CONTINUE 

56 


IPASS=IPASS+1 

57 

C 


60 

C 


61 

C 


62 


CALL. MXV<D2R.L1,L1) 

63 


CALL MXVCD2R>L2,L2) 
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CALL MXV<D2R.D1.D1> 

CALL MXV<D2R.XO. XO) 

C 

CALL CR0SS<D1 ,L1> 2U.2M) 

CALL CR0SS<2U,L1. YU>YH) 

XM=SQRT<L1<1)*>«‘2+Ll<2>t.*2+Ll<3)**2> 

irCXM. LE.O . ) GO TO 2000 

DO 5000 1=1,3 

XU<I>=Ll<n/XM 

YU< I )=YU<: I )/YM 

2U<I)=2U<I>/2M 

DC<3*I-2>=XU<I) 

DC<3+I-l> = YU<n 
DC<3*I)=2U<I) 

5000 CONTINUE 

CALL MXV<DC, L2, TEMP) 

A=TEMP<2) 

B=XM 

C=TEMP< 1 ) 

S=A+B/'2. 

IF<S.LE.O . ) GO TO 2000 
UIT<1 )=S*RHO*H 
C Wl.lO 

C CENTROID OF TRIANGLE 
C 

XBAR< 1 )=<B+C)/3 . 

XBARC2>=A,"3 , 

XBAR<3) = 0 . 

C iJll . 11 

C INERTIA MATRIX IN PANEL COORDINATES 
C 

. TMP=UT<1)/18. 

PMIP< 1 JsTMP+A^+S 

PMIP<5) = TMP:ti<B*H«2-B*C+C**2) 

PMIP<9)=PMIP< 1 )+PMIP<5) 

PMIP<2) = -TMP*A>t'<B/2 . -C) 

PMIP<4)=PhIP<2> 

PMIPC3)=0. 

PMIP<6)=0 . 

PMIP<7)=0. 

PMIP<8)=0. 

C Un . 12 

C ROTATE INEP-^TA MATRIX TO PARALLEL TO REFERENCE FRAME 

C 

CALL MXMSXSCPMIP, DC, TEMP) 

CALL XPOSE<DC, TEMP2) 

CALL MXM3X3<TEMP2, TEMP,PMI ) 

C UI1.13 

C COMPUTE CENTROID POSITION IN REFERENCE FRAME 

C 

CALL MXV<TEMP2,XBAR,TEMP) 

XBRFC 1 )=TEMP< 1 )+X0 < 1 ) 

XBRF<2)=TEMP<2)+X0 <2) 

XBRF<3)=TEMP<3)+X0 <3) 
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45 
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C UI1.14 

C VOLUME OF PANEL URT X-Y PLANE AND AREA ON Y-Z PLANE 
C 

V0L=0C<9)*S*XBRF<3) 

AREA=S'J‘0C<3> 

IFCAREA.LT.O . ) AREA=0 . 

C W1.15 

C MOMENT OF PANEL IN REF COORD 
C 

DO 5100 N=1 ,3 
UIT<N+l )=XBRF<N)*UJTa ) 

5100 CONTINUE 
C Itll . 16 

C COMPUTE MOMENTS AND PRODUCTS OF INERTIA IN REFERENCE FRAME 
C 

IXX=PMI < 1 > 

IYY=PMI<5) 

IZ2=PMI<9> 

IXY=-PMI (2) 

IX2=-PMI 

IYZ=-PMI<6> 

CALL PArCIXX.iUT<n,XBRF>UT<5)> 

C (lii . 17 

C SUBTOTAL OF MASS PROPERTIES 
C 

CALL TOTCSIXX, lt)T, 15) 

C W1 ■ 18 

C IS QUADRALATERAL COMPLETED ? 

C 

IF<IPASS.EQ. 1) GO TO 5020 
2000 CONTINUE 
- GO TO 80 
400 CONTINUE 

UJRITE< 10 . 6941 ) 

6941 FORMAK’ MADE IT TO 400’) 

C W1 . 19 

C SET DELTA CG ' S FOR NEUJ BASE LINE 
C 

1F<IBASE.EQ.0) GO TO 4041 
DO 4040 1=1,3 

4040 DCGCNO, I )=0 . 

IBASE=0 

C W1.20 

C START SUMMATION OF BLACK BOXES 

W 

4041 IRF=1 

403 CALL LOADCSIXX, TEMP, 15) 

409 IF<TEMP<1) .LT, 1 .E-4) GO TO 402 
DO 405 1=2,4 
J=14+I 

IF^ ASS<TEMP< I > ) , LT . 1 , ) TEMP<I)=0. 

TEriP<J) = TEMP'CI) 

405 TEMP< I )=TEMP( I )/TEMP< 1 ) 

CALL PAT<TEMP<5), TEMP<1),TEMP<2)/TEMP<5)) 
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67 UMP».5;=TCMP<5>-2i«<TEMP< i7)>*TCMP<3>+TEMP<18)*TEMP< 

70 TEIlP<6> = TEMP<6>-2>t<<TEMP< 16)+TEMPC2)+TEMP< 18>*TEMP< 

71 TEHP<7> = TEMP<7)-2 + <TEMP<16) + TEMP':2)+TEMP<17> + TEMP< 

72 TEMPC85 = TEMP<8?-TEMP< 16)><‘TEMP<3>-TEMP< 17>*TEMP<2) 

73 TEMP<9)=TEMP<9)-TEMPa6>*TEMP<4)-TEMP< 18)*TEMPC2) 

74 TEMP< 10 )=TEMP< 10 >-TEMP< 17)*TEMP<4)-TEMP< 18>*TEMP<3 

75. 402 IFUREFL.EQ.O .OR. IRF.GT.l) CO TO 401 

76 TEMP<3>=-TEMP<3) 

77 TEMP<8)=“TEMP<8> 

100 TEMP<10)=“TEMP<10> 

101 CALL PAT<TEMP<5>>TEMPa>,TEMP<2)/TEMP<5)) 


PAGE 


4)> 

4)> 

3)) 


) 


WftBCOM 


VER I 


REV A 


16 JUN 76 


PAGE 13 


1 



TEMP<2)=TEMPC2)*TEMP< 1) 

2 



TEHP<3)=TEMP<3>*TEMP<1 > 

3 



TEMP<4) = TEMP<4)H.TEMPa) 

4 



CALL TOT<SIXX,TEMP, 15> 

5 



IRF = 2 

6 



GO TO 403 

7 


401 

CONTINUE 

10 



IF<<BLKBX) .OR. (WTF.EQ.O . > .OR. <IRF.EQ.3>> 

11 

C 

UJl .21 

12 

C 

DENSITY CALCULATION FOR GEOMETRY OPTION 

13 

C 



14 



RHO=WTF/<TEMP< 1 1)*H> 

1? 



SIXXCl > = TEMPCn = WTF 

16 



DO 470 1=5. 10 

17 



SIXX< I) = TEHP< I ) = TEMP<n*RHO 

20 


470 

CONTINUE 

21 


469 

CONTINUE 

22 



DO 404 1=5. 10 

23 


404 

TEMP< I ) = TEMP< I )/4633 .056 

24 



FAC2=FACL*FACL 

25 



TEMP < 1 1 ) =TEMP C 1 1 ) /FAC2 

26 



TEM?a3)=TEMP0 3)/FAC2 

27 



TEMP< 12)=TEMP< 17)/(FAC2»FACL> 

30 



IF<IWRITE.EQ.0> GO TO 4999 

31 

0 

UJl .22 

32 

C 

WRITE RESULTS TO DISK 

33 

C 



34 



WRITEC20 . 1000 > NMB < 1 > . NMB < 2 > . <TEMP<I>. 1=1 

' 35 



WRITE<20 . 10 1 1 ) <TEMP<I). 1=11,135 

36 

4999 

IFC<NFLAG<65 .EQ. 1 ? .QR. UFLP.EQ. 155 GO TQ 

37 



IFC<NFLAG<55 .EQ. 15 .AND. <IRF.EQ.355 GO TO 

40 

C 

W1 , 23 

41 

C 

WRITE RESULTS TO ANK OR PPL 

42 

e 



43 



WRITECNUMO. 1000 5 NHB< 1 5 . NMB<25 , CTEMP<I5, I 

4 4 



WRITE<NUMO. 10115 <TEMPCI>, 1=11, 135 

45 



IFCNUMO.EQ. 135 CL0SE<NUM05 

46 



IF< IRF.EQ.35 GO TO 407 

47 



NMBBC 1 5=NMB< 1 5 

50 



NMBB<25=NMB<25 

51 


407 

IF<<ELKBX5 .OR. CIRF.EQ.355 GO TO 410 

52 



CALL LOAD<TEMP. WT, 155 

53 


410 

IF<IRF.NE.35 GO TO 411 

54 



WRITE<10, 69915 NFLAG<65 

55 

6991 

FORMATC NFLAGC65=‘ , 13) 

56 



IF<NFLAG<65 . EQ . 15 RETURN 

57 



IF< IFLP . EQ . 1 5 GO TO 360 

60 



PAUSE 

61 



RETURN 

52 


41 1 

IF : . NOT . SUB5 GO TO 333 

63 



DO 334 JJ=1,13 

64 

C 

U:1 

.24 

65 

C 

SUBTRACT COMPONENT'S MASS PROP. FROM TOTALS 

66 

C 

- 



GO TO 469 


410 

410 


88 


uf^ecoM 


VER 1 


REV ft 


16 JUN 76 


PACE 13 


67 

70 

71 

72 

73 

74 

75 

76 

77 
100 
101 
102 

103 

104 

105 

106 
107 
110 
111 
112 

113 

114 

115 

116 
117 
120 
121 
122 

123 

124 

125 

126 
127 

130 

131 

132 

133 

134 

135 

136 


334 C J T)=-SIXX< JJ) 

333 CONTINUE 

CAUL TOKTIXX. SIXX, 15> 
IF<NFLftC<4) .NE, 1) GO TO 335 
NAME C 1 > NO>=NMe< 1 ) 
NAME<2,N0)=NMB<2) 

335 IFCIFLP.EQ.U 60 TO 336 
ROOT=*STNAM£ 

ENCODERS, 1700 . NMB<0) ROOT 
ENCODECS. 1710 . NMP<2>) NTIME 

336 IRF=3 

CALL LOADCTIXX. TEMP. 15) 

GO TO 409 
C ID1.25 

C BLACKBOX CONTRIBUTIONS 
C 

’000 CONTINUE 

CALL LOADCUJT. SIXX. 15) 

DO 70 ir 1=1.10 
FAC=FftCI*4633 .056 
IFCI .LT.5) FAC=FACG 
IF<I .EQ. 1) FftC=FACW 
SIXX< I ) = SIXXd )*«FAC 
7010 CONTINUE 

SIXX<2)=SIXX<2)+DELX( 1 ) 

SIXXC3)=SIXX<3)+DELY 

SIXX<4)=SIXX<4>+DEL2 

CALL MXV<D2R. SIXX<2) .SIXXCa) ) 

PMIPC 1)=SIXX<5) 

PMIP<5)=SIXX<6) 

PMIP<9)=SIXX<7) 

rMlP<2)=-SIXX<S> 

PMIP<4>=-SIXX<8) 

PMIP<3)=-SIXX<9) 

PMIP<7)=-SIXX<9) 

PMIPC6)=-SIXXC 10 ) 

PMIP<8)=-SIXX< 10 ) 

CALL XP0SE<D2R. TEMP) 

CALL MXM3X3<PMIP. TEMP. TEMP2) 
CALL. MXM3X3<D2R. TEMP2.PMI) 
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SIXXC5)=PMI<1) 

SIXX<6>=PMI 
SIXX<7?=PMI C9) 

SIXX<8)=-PMI<2) 

SIXX<9)*-PMI<3) 

SIXXC10>=-PMI<6) 

CALL PAT<SIXX<5). SIXX<1>,SIXX<2>.SIXX<5>> 


10 


SIXX<2)=SIXX<:2)>t‘SIXX< 1) 

11 


SIXX<3)=SIXXC3)*SIXX< 1 ) 

12 


SIXX<4)=SIXX<4)>4«SIXX< 1) 

13 


WRITEC 10 , 6900 ) 

14 


6900 FORMAK ' GO TO 400 ' > 

15 


GO TO 400 

16 


IMAGE GEO 

17 

C 

QJ2.0 

20 

C 

DISPLAY GEOMETRY INPUTS 

21 

C 


22 


CLEARPEN<1 > 

23 


CHAR <”9B9X90 559YSi0 20 GEOMETRY 

24 


CHAR<G1N) 

25 


PEN0N<1) 

26 


MENUC "9B9X0O2O8Y9O3ONAME 


INPUTS") 


=9 . 8A9 
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1 * 8X8020SYSCI34INPUT PVV »8,8A9 

2 * 8X80208Y9040GEOM FILE =8.9A8 

3 ♦ 9X8020&Y9044RHO =8.9A9 

4 * SXS0208Y9050H =8.9A9 

5 ♦ 8X90209Y9054FACL -8.9A9 

6 * 8X9U209Y9060FACI sq . si^9 

7 * 8XS0208Y9064FACU »8.9A9 

10 ' * 9X90209Y9070FACG =9.9A8 

11 * 8X90208Y9074FIX UT =9 . 9A9 

12 * 9X81 10SY9030 IREFL =8,8A8 

13 * 8X81 108Y9034IIURITE =8.8A8 

14 ♦ 8X81 108Y9040SUB =8 . 8A8 

15 * 8X81 108Y9fJ44C>ELX =8 . 8A9 

16 ♦ 8X91 108Y9050DELY =8 . 9A9 

17 4. 9X81 108Y9054DEL2 *9 , 9A8 

* 8X81 108Y9060PSI =9 . 9A9 

21 * 9X91109Y8064THETA =9 , 8A9 

22: » 8X81 108Y8070PHI =9 . 9A9 

23 * 8X81 108Y8074BASE CHNGE=8 . 9A9 

24 * 8X80709Y9104PROCEED9 . 92" ) 

25 PEN0FF<1> 

26 GO-TO <300,30 1,302,303T-304^305i306i^07v308/298,309/ ■ 

27 * 310.311.312.313>314.315>316/317,297i299> 

30 GO TO 500 

31 300 IFLAGC3)=. ,RUE. 

32 IFLAGCSl )=.TRUE. 

33 IFLAGC30)=.TRUE. 

34 GO TO 500 

*35 301 IFLAG<4)= . TRUE . 

36 IFLAGCSO )= , TRUE . 

37 GO TO 500 

40 302 1FLAG<5)=. TRUE. 

41 ’ IFLAG<31 )=. TRUE. 

42 IFLAGC30>=.TRUE, 

43 GO TO 500 

44 303 IFLAG<6>=, TRUE. 

45 IFLAG<30)=.TRUE. 

46 GO TO 500 

47 304 IFLAG<7)=.TRUE. 

50 IFLAG<30)=.TRUE. 

51 GO TO 500 

52 305 IFLAG<S>= . TRUE . 

53 IFLAe<30 )= , TRUE . 

54 GO TO 500 

55 30 6 IFLAG<9>fa: . TRUE . 

56 IFLAG<30)= . TRUE . 

57 GO TO 500 

60 307 IFLAG<10)=.TRUE. 

61 IFLAG«C30> = .TRUE. 

62 GO TO 500 

63 308 IFLAG< 1 1 >=.TRUE, 

64 IFLAG<30)=.TRUE. 

65 GO TO 500 

66 309 IFLAG< 12)=.TRUE. 
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67 

70 

71 

72 

73 

74 

75 

76 

77 
100 
101 
102 

103 

104 

105 

106 
107 
110 
111 
112 

113 

114 

115 

116 
117 
120 
121 
122 

123 

124 

125 

126 
127 

130 

131 

132 

133 

134 

135 

136 

137 

140 

141 


C 

C 

C 


irLAG(30)=.TRUE. 

CO TO 500 

310 IFLAG<13>=.TRUE. 
IFLAGC30)= . TRUE. 

CO TO 500 

311 IFLAG<14)=.TRUE. 
IFLAG<31)=.TRUE. 
IFLAG<30>=.TRUE. 

GO TO 500 

312 IFLAGa5)=.TRUE. 
IFLAG<30)=. TRUE. 

GO TO 500 

313 IFLAG<16)=.TRUE. 
IFLAG<30>=. TRUE. 

CO TO 500 

314 IFLAC<17>=.TRUE. 
IFLAGCSOJs.TRUE. 

GO TO 50 0 

315 IFLAG<1S)=.TRUE. 
IFLAG<30)=.TRUE. 

GO TO 500 

316 IFLAG<19>=.TRUE. 

I FLAG <30 TRUE. 

GO TO 500 

317 IFLAG<20)=.TRUE. 
IFLAG<30)=,TRUE. 

GO TO 500 

298 IFLfiG<22>= , TRUE. 
IFLAGOO >= . TRUE. 

GO TO 500 

297 IFUAG<23>= , TRUE . 

. IFLAG<30>= , TRUE. 

GO TO 500 

299 IFLAG<21 )= . TRUE . 
IFLAG<30>=,TRUE. 

500 CALLS I MG GRIDW 
RETURN 
IMAGE BKBX 
Ul,' . 0 

DISPLAY BLACK BOX INPUTS 


CLEARPEN<1 > 

CHAR< "aBaXaOBOaVaOBOBLACK BOX INPUTS") 
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1 CHAR<BIN> 

2 PEN0N<1) 

3 MENUC'aBQXaO 10&Y8030MAME= 9 . SAS 

4 ♦ 8XO0109YO034WT = 9.9A9 

5 * 8X90 109YO040XCG = 8.8A8 

6 <*• 8X90709Y9040YCG = 8 . 9 As 

7 ♦ 8X91509Y90402CG = 8.8A8 

10 ' * 8X80 108Y9044IXX = 8 . 8A8 

11 «*« 8X80708Y8044IYY =» 9 . 8A8 

12 * 8X81508Y9044IZZ = 9 . 9A8 

13 * 9X80108Y9050IXY = 9.8A8 

14 * 8X80709Y8050 lYZ = 9 . 9A8 

15 * 9X81508Y9050 IXZ = 9 . 9A9 

16 * 8X80 108Y8054ASF = 9 . 9A8 

17 ♦ 8X80709Y9054VOL = 8 . 9A8 

20 * 8X81509Y9054AFR = 8.9A8 

21 ♦ 9X90209Y9064FACL * 8.8A8 

22 * 8X90209Y9070FACI = 8 . 8A9 

23 f« 8X90208Y90 74FACW = 8.9A8 

24 ♦ 8X90208Y8100FACG = 9.8A8 

25 * 8X80209Y91040VL = 9.9A9 

26 ♦ 9X80208Y9H0 UIRITE = 9.9A8- 

27 * 8X80209YO1 14IPRINT = 8 . 8A9 

30 >♦« 8X90208Y9120GEOM FILE = 9.9A9 

31 ★ 8X81 109V9064DELX = 9 . 9A8 

.*■ 32 * 9X9U08Y9070DELY = 8 . 9A8 

33 * 8X91 108Y9074DEL2 = 8.9A9 

34 ♦ 8X91 109Y8100PSI = 8 . 8A8 

' 35 * 8X91108Y8104THETA = 8.9A8 

36 * 9X81 108Y81 lOPHI = 8 . 8A8 

37 * 8X91 108Y91 14BASE CHNGE= 8.8A8 

40 f 9X80708Y9130PROCEED 8.82") 

41 PENOFFCl) 

42 60 TO <364> 365. 366> 367. 368> 369, 370 , 371 , 372- 373^ 374, 375> 376, 377, 

43 * 378,379,380,381,332,383,391,397,384,385,386,387,388,389,392,390) 

44 GO TO 393 

45 364 IFLAGC29)= . TRUE , 

46 IFLAG<30)= . TRUE. 

47 IFLAG<31)=.TRUE. 

50 GO TO 393 

51 ^ 365 IFLA6<3)=. TRUE. 

52 ' IFLAG<30)=.TRUE. 

53 GO TO 393 

54 366 IFLAG<4)=.TRUE. 

, 55 1FLA6<30)=.TRUE. 

56 60 TO 393 

57 367 IFI.AG<5)= . TRUE. 

60 IFLAG<30)=.TRUE. 

61 60 TO 393 

62 368 IFLAG<6)= . TRUE . 

63 IFLAG<30)=.TRUE. 

64 GO TO 393 

65 369 IFLA6<7)=. TRUE. 

66 - IFLAGC30)=. TRUE. 

93 




ItlABCUM 

VER * REV A 16 

67 ' 

CO TO 393 

70 

370 IFLAG<8)=. TRUE. 

71 

IFLAG<30)=.TRUE. 

72 

GO TO 393 

73 

371 IFLAG<9)=.TRUE. 

74 

IFLAG<30>=.TRUE. 

75 

GO TO 393 

76 • 

372 IFLAGC10)=.TRUE. 

77 

IFLAG<30>=.TRUE. 

100 

GO TO 393 

101 

373 IFLAG< 1 1 >= . TRUE . 

102 

IFLAG<30>=.TRUE. 

103 

GC TO 393 

10 4 

374 IFLAG<12>=.TRUE. 

105 

IFLAG<30)=.TRUE. 

106 

GO TO 393 

107 

375 IFLAG<13)=.TRUE. 

110 

IFLAGC30)=. TRUE. 

111 

GO TO 393 
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376 IFLAG<14>=.TRUE. 
IFLAC<30)=.TRUE. 
r<0 TO 393 

377 IFLAG<15)=.TRUE. 
1FLAC<30>=.TRUE. 
GO TO 393 

378 IFLAGC16>«=.TRUE. 
IFLAGC30 )= . TRUE. 
GO TO 393 

379 IFLAGa7)=.TRUE. 
IFLAG<30)=.TRUE. 
CO TO 393 

380 IFLAG<18)=.TRUE. 
IFLAG<30>=.TRUE. 
GO TO 393 

381 IFLAG<19)=.TRUE. 
IFLAG<30)=. TRUE. 
CO TO 393 

382 IFLAG<20>=.TRUE. 
IFLAGC30>=. TRUE. 
IFLAG<31>=.TRUE. 
GO TO 393 

383 IFLAG<21)=.TRUE. 
IFLAGOO >=. TRUE. 
GO TO 393 

384 IFLA6C22?=. TRUE. 
IFLAGOO >=.TRUE. 
GO TO 393 

385 IFLAG<23)=.TRUE. 
IFLAG<30)=.TRUE. 
GO TO 393 

386 1FI.AG<24) = .TRUE. 
IFLAG<30>=.TRUE. 
GO TO 393 

387 1FLA6(25)= . TRUE. 
IFLAG<30)=.TRUE. 
GO TO 393 

388 IFLAG<26)=.TRUE. 
IFLAG<30)=.TRUE. 
GO TO 393 

389 IFLAG<27>==.TRUE. 
IFLAG<30)=.TRUE. 
GO TO 393 

390 IFLAG<2S>=,TRUE. 
IFLAGOO >= . TRUE. 
GO TO 393 

391 IFLAG<32)=.TRUE. 
IFLAG<30)= . TRUE. 
GO TO 393 

392 IFLAG<33)=. TRUE. 
IFLAG<30>=.TRUE. 
GO TO 393 

397 IFLAG<34>=. TRUE. 
IFLAG<31)=.TRUE. 
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67 

70 

71 

72 

73 

74 

75 

76 

77 

too 

101 

102 

103 

164 

105 

106 
107 
110 
111 
112 

113 

114 

115 

116 
1 17 
120 
121 

-122 

123 

124 


IFLAG<30)«.TRUE. 

393 CALLSIMG GRIDUl 
RETURN 

SUBROUTINE RESET 
C UI4.0 

C INITIALIZE SUBROUTINE VARIABLES AFTER EACH INPUT 
C ; 

IINUM»0 
RRNUM=0 . 

DO 395 1=1.40 

395 IFLAGC I >=. FALSE. 

RETURN 

SUBROUTINE RE3ET69 
C U5.0 

C INITIALIZE SUBROUTINE VARIABLES AFTER EACH CALCULATION 
C 

IF<NFLAG<5) .NE. 1) GO TO 396 

NM8<1 )=5HDATA 

NMB<2>=5HFILE 

IDPV=110 

IPV=110B 

IFILE<n=5H 

IFILE<2>=5H 

396 IREFL=HURITE=IPRINT=0 

WtF = DELX=DELY=DELZ=PHI=THETA=PSI=0 . 
RHO=H=FACI=FACU)=FACG=l . 

FACL=12. 

SUe=. FALSE. 

RETURN 

END 
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SUBROUTINE IDABLST 
C Ll.O 

C ; MENU COMPUTATION FOR COMPONENT DISPLAY 
C ' 

COMMON/'PARA/' LFLG< 10 > / MFLAG<4) / SUM / LETTER. NO, NFLG. INO. NFLAG < 12) 

•, / JJJ/ KKK/ LLL/ MMM/NNN, NAME<2/ 60 ) / NME<2> > I INUM, RRNUM 
COMMON/'CNT/ IPASS. NUP. NAMER<3> / IPV. ISTRT, lEND. IT. LINES * 
COMMON^OPT/ NHIT/OPDPL 

COHMON/STR/ NTOT / VALUE 1 < 60 > / VALUE2C60 > / VALUES C 60 ) / VALUE4<60 ) / 

♦ VALUES <60 ) / VALUE6<60 ) / VALUE7C60 ) /VALUES <60 ) / VALUES <60 ) / 

*; VALUE10<60) / VALUEl 1 <60 ) / VALUE12<60 ) / VALUE13<60 ) 

COMMON/liiB/' NMB<2)/WT<15)/TEMP<18)/ lUJRITE/ IJ. ITNM<2) / TIXX< 15) / NUMO 

* / IPRINT/PI»»T<3)/DCG<60/3)/ IGNME<3/60) 

DIMENSION TOTK20)/TOT2<20)/TOT3<20>/TOT4r20)/TOT5<20) 

INTEGER MNU 1 < 20 0 ) / MNU2 < 200 )/ MNU3 <200)/ MNU4 < 200)/ MNU5 < 20 0 ) 

INTEGER FMT<8) ,MASK<3) / NAME/DX/DY,iPUJD4/NME/MNU6<200) 

REAL VALUEl/ VALUE2/VALUE3/VALUE4 
DATA HF/10/' 

DATA MASK/'8*1H /' 

C 12345678S0 1234567890123456789012345678901234567890 

DATA FM!/' <9H<"9SSiBaXa/03,3H8V3/03/ / 6HQ . SA " ) ) ' z' 

IF<NUP .EQ. 1 ) GO TO 100 
IF<IPASS. NE.O) GO TO 900 
• 10 IPr.S3=lPASS+l 
C L* . 1 

C COMPONENT'S NUMBER MENU 
C 

DY=24 

DX-6 

LABEL <MNU6) 

2SET<0.0F) 

DO 200 J=ISTRT/IEND 
' FMT<8) = 5H9. ")) 

IF<J.EQ, lEND) FMT<8)=5H3A">) 

FMT<6)=3H1I3 

ENCODE<40/FMT/MASK/H) DX/DY/J 
DY-DY+4 

URITE<16/MASK) 

200 CONTINUE 
C LI .2 

C COMPONENT'S NAME MENU 
C 

DY=24 

DX=10 

LABEL<MNU1) 

ZSET<0,0F) 

. DO; 300 J=ISTRT/ lEND 
FMt<S)=5HaA”)) 

IFtJ. EQ. lEND) FMT<8)=5Ha2")) , 

FMT<6)=3H2A5 

ENrODE<40/ FMT/MASK/H) DX/ DY. NAME< 1 / J) / NAME< 2/ J) 

DY * DY+4 
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i 

i 

4 

s 

i 

7 

to 

11 

12 

13 

14 

15 
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33 
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44 

45 
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62 
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65 

66 


16«MftSK> 

300 corniNUC 
F.M>CIST 
C L» .3 
C TOTALS- MENU 
C 

DY*DY+4 
LABELCTOTl > 

ZSETCO .0F> 

ENCODE<40,FMT/MASK.H) DX.OY/ ITNM< 1 ) / ITNMC2) 
mRITEU6.MASK) 

ENDLIST 
C Ll-4 

C COMPONENT’S WEIGHT MENU 
C 

DY«s24 

DX«25 

LAjjEL<MNU2) 

2SET<0.0F> 

00 400 J=ISTRT. lEND 
FMT<8)=5HaA">> 

IF<J.EQ . IEND> FMT<8)=5H8Z">) 

FMT<6>==5H 

CALL DECFMT<VALUEKJ>,NF.FMT<6>) 

ENC0DEC40> FMT>MftSK/H) OX. DY> VALUEl < j) 

DY=DY+4 

WRITE<16>MASK) 

400 CONTINUE 
ENDLIST 
C LI. 5 

C T:TAL WEIGHT MENU 
C - 

' DY=DY+4 

LABELCT0T2> 

2SET<0 .OF) 

ENCODE C 40, FMT> MASK. H) OX. DY. TEMP< 1 > 
WRITE<16,MASK) 

ENDLIST 
C LI. 6 

C COMPONENT'S X-CG MENU 
C 

DY=24 

DX=40 

LABEL <MNU3) 

2SET<0 ,0F) 

DO 500 J=ISTRT . lEND 
I FMT<3)=5H8A" ) ) 

I IF<J.EQ. TEND) FMT<8)=5Hq2" ) > 

i FMT<6>=5H 

! CALL DECFMT<VALUE2<J).NF.FMT<6)) 

ENCODE<40,FMT.MASK.H> DX. DY, VALUE2C J> 

DY=DY+4 

WRITECS.MASK) 

500 CONTINUE 
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124 
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ENDLIST 
C : LI .7 

C ; TOTAL X-CG MENU 
C 

0Y=DY+4 
LABEL CT0T3) 

•2SET<0 .OF) 

ENCODE<40.FMT.MASK>H) DX. DY> TEMP<2> 
UIRITEC16.MASK) 

ENDLIST 
C wl.8 

C COMPONENT'S Y-CG MENU 
C 

0Y=24 

DX=55 

LABEL<MNU4) 

ZSET<O.OF) 

DO 600 J=ISTRT. lEND 
FMT<8)=5HaA")> 

IFCJ.EQ. lEND) FMT<8)=5H92">) 

FMT<6>=5H 

CALL DECFMT<VALUE3<J) »NF/FMT<6>) 
ENCODE < 40. FMT. MASK. H) DX. DY. VALJE3< J) 
DY=DY+4 
WRITE<1£.MASK> 

600 CONTINUE 
ENDLIST 
C LI. 9 

C TOTAL Y-CG MENU 
C 

DY=DY+4 

LABEL<T0T4) 

* 2SET<0 .0F> 

ENCODE < 40. FMT.M ASK. H) DX. DY, TEMPC3> 
WRITE<16.MASK) 

ENDLIST 
C Ll.lO 

C COMPONENT'S 2-CG MENU 
C j 

DYi24 

DX=70 

LABEL<MNU5> 

ZSETCO.OF) 
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to 700 J=ISTRT,IEND 
rHT<8)=5HSA")) 

IFCJ.EQ. lENO) FMT<8)=5HSZ")> 

Ff1T<6>=5H 

CALL DECFMT<VALUE4<J)»HF.FMT<6)) 
ENCODE<40/FMTiMASK.H) DX, DY. VALUE4< J> 
DY=DY+4 

UIRITE<16,MASK> 

700 CONTINUE 
ENDLIST 
C Ll.ll 

C TOTAL Z-CG MENU 
C 

DY«DY+4 

LABEL<T0T5> 

ZSET<0.0F) 

ENCOOE<:40 / FMT.MASKiH) DX. DY, TEMP<4) 

UIRITEC 16, MASK) 

ENULIST 

NUP=0 

800 DISPLAY UABMHU 
900 1F<LFLG<3) .NE. 1) RETURN 
IFXJJJ.EQ.O) GO TO 1015 
NAMECl , JJJ)=NME< 1) 

NAME<2/ JJJ)=NME<2) 

NUP=1 

RETURN 

1015 IF<KKK.EQ.O) GO TO 1030 
VALUEKKKK) = RRNUM 
NUP=1 
RETURN 

1030 IF'LLL.EQ.O) GO TO 1045 
P»UT<1)*VALUE2<LLL) 

VALUE2<LLL)=RRNUM 

DCG<LLL> 1 )=DCG<LLL, 1 ) + <VALUE2<LLL)-PIUT< 1 ) > 

NUP=1 

RETURN 

1045 1F<MMM.EG!.0) GO TO 1060 
PUIT<2)=VALUE3<MMM) 

VALUE3<MMM^=RRNUM 

DCG<MMM,2)=DCG<MMM,2) + CVALUE3<MMM)-PUIT<2)> 

NUP=1 

RETURN 

1060 IFCNNN.EQ.O) GO TO 800 
PUIT(3)=VALUE4<:NNN) 


100 


WABLST 


VER 1 


REV ft 


PAGE 4 


16 JUN 76 


1 

2 

3 

4 

5 

6 
7 
10 
11 
12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 
/ 32 

33 

34 
* 35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 

. 55 

56 

57 
60 
61 
62 

63 

64 

65 

66 


VALUE4(NNN)=RRNUM 

DCG<NNN>3)-DCG(NNN«3>*^<VALUE4(NNN)-PUIT(3>> 

NUP=*1 

RETURN 

IMAGE (UAeMNU 
C L2 . 0 

C DISPLAY COMPONENTS’ MASS PROPERTIES 
C 

CHAR<"SS0XSO1OOYSO2O COMP WT X-CC 

* Y-CG Z-CC") 

CHAR<MNU6) 

CLEARPEN<1) 

♦PWD4=0 

PENONCl) 

MENU<MNUO 

PEN0FF<1) 

CHAR<T0T1> 

IF<<LFLG<10) .EQ. 1) . OR. <HFLAG<1) ,EQ. n .OR. <MFLAGC2) ,EQ. O 
f .OR. <NHIT. EQ. 1 >) GO TO 1101 
IF<4PWD4. EQ. 0 ) GO TO 1101 
LFLG<6)=1 

JJJsNO=tPUJD4.*<-IT.>t'LINES^ .. 

MFLAG<1)=1 
1101 CLEARPEN<1> 

♦ PUJD4 = 0 
PENONCl) 

MENU<MNU2> 

PEN0FFC1> 

CHAR<T0T2> 

IFC<LFLG< 10) ,EQ. 1 ) . OR . <MFLAG<1) .EQ. I ) . OR . <MFLAG<2) . EQ . 1) 

* .OR..<NHIT.£Q.D) GO TO llll 
IF<iPlUD4.EG!.0) GO TO 111! 

LFLG<6)=1 

KKK=N0='$PtUD4 + < IT’^'LINES) 

IF<KKK.NE.O) IN0=2 
MFLAG<2)=1 
I: ii CLEARPEN<1) 

♦ PUJD4=0 
PEN0N<1) 

MENU<MNU3) 

PEN0FF<1> 

CHAR<T0T3) 

IF<<LFLGC10) .EQ . 1) .OR. <MFLAG<1> .EQ . 1 ) . OR . <MFLAG<2) . EQ . 1 > 

* .OR. <NHIT.EQ. 1)) GO TO 1121 
IF<fPWD4.EQ.O) GO TO 1121 
LFLG<6)=1 

LLL=N0=-rPWD4+< IT^LINES) 

IF<LLL.NE.O) IN0=3 
MFLAG<2)=1 
1121 CLEARPENO) 
iPW04=0 
PENONCl) 

MEHU<MNU4) 

* PENOFFCl) 


101 
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67 

70 

71 
7g 

73 

74 

75 

76 

77 
100 
101 
102 

103 

104 

105 
1C6 
107 
110 
111 
112 

113 

114 

115 

116 
117 


CHARCT0T4> 

rF<<LFLCU0> .EO. D.OR. <MFLAG<1) .EQ.1>.0R.<MFLAC<2>.EQ. 1> 

• .OR. <NHIT.EQ. 1 )> CO TO 1131 
lF<t?aiD4.EQ.O) GO TO 1131 
i.FLG<6> = l 

MMM*N0=tPU)04+< IT*LINES> 

IFCMMM.NE.O) IN0*4 
MFLAG<2)*1 
1131 CLEARPENd) 

♦ PUID4=0 
PEN0N<1) 

MENU<MNU5> 

PENOFFCl) 

CHAR<T0T5) 

IF<<LFLG<10) .EQ, 1 ) . OR . CMFLAG < 1 > . EQ . I > . OR . <MFLAG<2> . EQ . 1 > 

* .OR. <NHIT.EQ. 1)) GO TO 1200 
IF<tPWD4.E(3.0> GO TO 1200 
LFLG<6>=1 

NNN=NO=f PUID4+< IT^LINES) 

IFCHNN.NE.O) IN0=5 
MFLAG<2)=1 

1200 IF<<MFLAG<1> .NE.O) .OR. <MFLAG<2> .NE.0>) NFLG»1 
CALLSIMG OPTIONS 
RETURN 
END 
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1 

2| 

3 

4 

5 
6; 
7 
10 

12 

13 

1?^ 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 


SUBROUTINE UABSET 
C Tl.O 

C ARRAY SETUP FOR UABCOM 
C 

COMMON^STR/ NTOT, VALUEK60 ) , VALUE2<60 > * VALUE3<60 ) . VALUE4<60 ) , 
f VALUE5<60 > , VALUE6<60 ) , VALUE7<60 ) , VALUE8<60‘) , VALUE9<60 > , 

♦ VALUE10<60) , VALUElt <60), VALUE12<60>, VALUE13<60) 

COMMON/PARA/ LFLGdOS.MFLAGM), SUM, LETTER/NO, NFLG. IN0/NFLAG<12>> • 

; <•« JJJ,KKK,LLL,MMM/NNN,MAME<2,60)/NME<2), IINUM/RRNUM 

COMMOH-^WB/' NMB<2?,UT<15),TEMP«:iB), lURITE, IJ, ITNM<2) , TIXX< 1 5) , NUMO 

* , IPRIHT,PUJT<3)/DCC<60,3>, 1GNME<3,60) 

DIMENSION V<60, 13> 

EQUIVALENCE < V <1 , I > , VALUEK I > > 

DO 50 1=1,15 
! TIXXCI)=0. 
bo CONTINUE 

DO 100 IJ=l,NTOT 
I NMB<1) = NAME<1, IJ> 

NMB<2> = NAME<2, IJ) 

DO 200 IK=1, 13 
UT<IK>=V<IJ, IK) 

200 CONTINUE 

CALL UABCOM 
1,00 CONTINUE 

! I^OURITE.EQ. 1) CLOSE<20) 

ITHM<1)=5HT0TAL 

ITNM<2)=5HS 

RETURN 

END 
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1 

2 C 

3 C 

4 C 

5 C 

7 C 

10' C 

11 C 

12 C 

13 C 

14 

is; 

16 


t7 

20 

21 

22 

23 

24 
25i 
26 
27 
30 


3l! 

32 

33 

34 

35 


SUBROUTINE XFORM<A. C.X> 

THIS SUBROUTINE COMPUTES THE TRANSFORMATION MATRIX FOR A 
YAU«P1TCH«R0LL ROTATION SEQUENCE. 

INPUTS! A s YAlil ANGLE IN DEGREES 

B ■ PITCH ANGLE IN DEGREES 
C « ROLL ANGLE IN DEGREES 

OUTPUT! X * 3X3 MATRIX OF TRANSFORMATION DIRECTION COSINES 

DIMENSION X<9) 

SA*SINCA/57. 295779) 

CA*C0S<A/57. 295779) 

SB=SIN<B/'57. 295779) 

CBsC0S<B/57. 295779) 

SC=SIN<C/57, 295779) 

CC=C0S<C/57. 295779) 

X<1)=CB>*<CA 

X<2)=SC*SB*CA--CC>!<SA 

X<3)=CC*SB*CA+SC+SA 

XC4)=CB+SA 

X<5)=SC*SB>!‘SA+CC>^CA 

X<6)=CC*SB*SA'SC*CA 

X<7)=-SB 

XCS)=SC*CE 

X<9)=CC*C‘j 

RETURN 

END 
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1 ! 

s! c 

t c 

4 ! C 

5 

6 
?! 

10 

11 

12 

13 

14 

15 

16 
17 
20 


SUBROUTINC XP0SC<A.B) 

6 * TRANSPOSE OF A 

DIMENSION A<9),8<9> 
BCn-A<l> 

B<2>=AC4> 

BC3)=A<7) 

B<4)=A<2> 

B<5>=A<5) 

B<|6>=A<8> 

B<|7)=A<3) 

B<|8)=A(6> 

B<9>»A<9) 

RETURN 

END 


A S. B ARE 3X3 MATRICES 


